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ABSTRACT 


A computer program is P«.«t.d 

U free space, ^ 

finite conductivity and lumped loa._- a nd gain. The 

current distribution, impedance, radiation eT.icienuy 
program uses sinusoidal bases and balerkin s nx hod. 
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I. INTRODUCTION 


Reference 1 presents a computer program and reference 2 
presents the theory for thin-wi re antennas and scatterers in a horo- 
genecus conducting medium. The present program differs from 
reference 1 only in the following details: 

a. The exterior medium is- free space. 

b. The antenna is situated over a perfect ground plane. 

c. The wires have no dielectric sleeves. 

d. The frequency is real. 

e. Scattering problems are not considered. 

To avoid unnecessary duplication, it is assumed the reader is 
familiar with the program in reference 1. 

The program handles antennas constructed of straight wire 
segments. One or more segments may connect to the ground plane, or 
the antenna may be situated some distance away from the around 
plane. No segment has both endpoints on the ground plane. This 
program can readily be modified to handle core general situations 
as in reference 1. The program uses the delta-gap model for the 
generators . 

The method of images is employed to reduce the problem to a 
thin-wire structure in free space. The user sets up the geometry 
of the real wire configurat ion, and the pronram automatically sets 
up the image. If unlimited storage were available, one might set up 
a large mutual-impedance matrix for the wire antenna and its image 
in free space. Instead, this program takes advantage of the ground- 
plane symmetry and sets up the compressed matrix C(I,J). Only the 
currents on the real segments are treated as independent unknown 
quantities, and the image currents are dependent. All the currents, 
however, are influenced by the mutual couplings among all the seg- 
ments including the images. In taking advantage of tne mirror 
symnetry, we lose the advantage of havinq a sy rune trie matrix. 
However, the net gain in computational speed and storage is 
substantial . 

In practice, many wire antennas operate over a ground plane 
with finite conductivity and finite extent. In many cases, however, 
one may substitute a perfectly conducting ground plane of infinite 
extent without unduly disturbing the antenna current distribution or 
impedance. After these quantities have been calculated, one may 
then take into account the finite ground plane in calculating tne 
efficiency and patterns. The present program, however, assures an 
idealized ground plane throughout. 

The remaining sections present the computer program with 
enough explanation to enable an experienced engineer to use it. 
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II. THE INPUT DATA AND SUBROUTINE IMRE 


Fig. 1 is a Fortran listing of subroutine I WIRE. This 
subroutine is used to set up the input data. The following data 
must be read or programmed in I WIRE : 

A l wire radius a/\ 

CMM wire conductivity in megamhos/m 

DPH increment, in far-field angle ; in degrees 

FMC frequency in MHz 

TH elevation angle c in degrees for far-field pattern 

To define the shape and size of the wire antenna* the input 
data includes the coordinates XC(I), YC ( I ) , and ZC ( I ) of the wire 
endpoints, terminals and other current sampling stations along the 
wire axis. The unit of length is selected by the? user, and SCALE 
is the conversion factor such that XC(I)*SCALE gives the coordinate 
of point I in meters. NPGP denotes the number of points on the 
ground plane, and NRP is the number of real points including those 
on the ground plane. Coordinates are supplied only for the real 
points. The ground plane coincides with the xy plane. 

NRS denotes the number of real segments* and NSGP is the 
number of real segments having an endpoint on the ground plane. 

For each real segment J, the input data specifies the endpoints 
I A ( J ) and I B ( J ) . In assigning index numbers to the segments* the 

lowest numbers must refer to those having an endpoin t on the ground 
plane. In assigning index numbers to the points, the lowest numbers 
must refer to those on the ground plane. 

Set IWRCJ = 1 to obtain a writeout of the antenna currents; 
otherwise IWRCJ = -1. Set IWRITE = 1 to obtain a writeout of the 
antenna geometry; otherwise IWRITE = -1. If IN T = 0, the rigorous 
closed-form expressions will be used for the mutual impedance of 
sinusoidal dipoles and the calculations will tend to be slow but 
accurate. If 1 NT is a positive integer, Simpson's rule will be 
used for the mutual- impedance calculations. The closed-form ex- 
pressions are always used automatically for the most critical 
impedances. We usually use INT = 0 for multi -turn loop antennas 
with closely-spaced turns, and INT = 4 for general purpose. 

Simpson's rule uses INT integration intervals. Thus the accuracy 
and the execution time tend to increase with larger values of INT. 

Set NLD - 0 if there are no lumped loads; otherwise NLD = 1. 
ZLD(J) denotes the impedance (in ohms) inserted in segment J at 
endpoint I A( J) . A lumped load at endpoint IB(J) is denoted by 
ZLD(J + \RS). The user sets up only the real generators and lumped 
loads, and the program takes care of the images. 
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Fig. 1. Subroutine I W I RE . 


3 


REPRODUCIBILITY OF THE 
ORIGINAL PAGE IS POOR 


The generator location is defined by JGN. The numbering 
system for the generators is the saae as for lumped loads. If the 
generator is to be inserted in segrrent J at endpoint JA(J), the 
generator index JGN is the same as the segment index J. To insert 
a generator at endpoint IB(J), set JGN = J 4 NRS . VG(JGN) denotes 
the complex voltage of the generator. The reference direction for 
these voltages is from IA(J) toward IB(J). If the antenna is fed 
with several generators, delete JGN and merely input the generator 
voltages VG. 

NP and NS denote the n untie r of points and segments, re- 
spectively, for the complete system (antenna and image) in free space. 


III. DIE MAIN COMPUTER PROGRAM 

The main computer program is listed in Fig. 2. This p r0 9 r<3n 
calls subroutine 1WJRE for the input data. Then it calls subroutine 
1S0RT to generate and store' the data for the image points and image 
segments and the length DC(J) of each segment. Then ISORT generates 
a list of sinusoidal dipole nodes for the complete system (antenna 
and mage) in free space. Dipole node I has segments JA(I) and 
JB(I), terminals at point 12(1), and endpoints 11(1) and 13(1). 

This subroutine also generates the following information: 

ND(J) number of dipole codes sharing segment J 
MD( J , K) list of dipoles sharing segment J 
NCM size of the compressed matrix 
N number of dipole modes on the complete system 

The following quantities must be specified in the main program: 

ICC dimension for the compressed matrix C(I,J) 

ICJ dimension related to number of dipole modes N 

INP dimension related to number of points NP 

INS dimension related to number of segments NS 

In Fig. 2, all quantities having the sari? dimensions are dimensioned 
in the san? or adjacent staten?nts. The numerical values assigned 
to ICC, ICJ, INP and INS must agree with the dimensions actually 
reserved for the correspond' ng quantities in the COMPLEX and 
DIMENSION statements. ICC, ICJ, INP and INS must be at least as 
large as NCM, N, NP and NS, respecti vely . In Fig. 2, the main 
program is dimensioned for up to 150 .modes, 90 points, 100 segments, 
and a compressed matrix as large as 30 by 30. If the wire antenna 
makes no contact with the ground plane, the compressed matrix will 
bo exactly half a^ large as the full matrix. Otherwise NCM is some- 
what larger than N/2. 
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X(I),Y(I) and Z(I) ienote kx, ky and kz for point I, where 
k - 2v/\. If calculations are desired for a given anter.na jt severe' 1 
frequencies, t h e frequency DC LOOP will begin just below the call to 
I SORT. 


The main program calls subroutine IDANT to generate the corn- 
pressed open-circuit impedance matrix C(i,J). Th^n si routine AMI 
is called to obtain tr.e current 4 i > r ributicn CJ(I) find t*e radiation 
efficiency IFF. ANTI also calculates the complex power nput to the 
c;uc<>na, denoted by Y*l, and the t me-average input power G. If the 
antenna has only one generator anu VG(JCN) = (l.,D.), mor. Yll and 
Zll denote the antenna admittance and impedance* respect: vely. 

Finally, the antenn.i pattern is obtained by calling subroutine 
IfftD. th and PH denote the spnerical coordinates ■ and ; (in 
degrvos) of the distant observer. GPP and G7T denote tne .‘-polarized 
and -polarized power gains, respecti vely, ana DBP ana D3T are che 
decibel versions. The user may want to increment as well as * , 
but tfiis will require only a trivial change in the rain proqram. 

IFFLD is called once for each look angle (•■*?). Wnen - 9C - , GpP 
will vanish if the program has set up a valid system of images. 


IV. AN EXAMPLE 


Fiq. 3 shows a single antenna and its image, with a cotted 
line to indicate the ground plane. In Fig. 1, subroutine IWIRE s^ts 
up trie following input data for this antenna: 


VG(4) = 1 
NL D = 0 
CMM = 1 . 
NSGP = 2 
NpGP = 2 
NRS = 4 
NR? = 5 


unit voltage generator at endpoint IA of segment 4 
no lumped load: 

the wire conductivity is 1 nega * o/m 
2 real segments connect to the ground plane 
2 points on the ground plane 

4 real segments 

5 real points 


This planar antenna has 8 points (NP = 8) and 8 secrionts (NS = 8), 
and fhe numbering system is shown in Fin. 3. Note tr.at the lowest 
numbers are assigned to the two points on the ground plane and the 
two segments terminating on the ground plane. Tre points and seg- 
ments must bo labeled with consecutive positive integers 1,2, 3, . 
For a given segment J, it rakes no differvn ce which end is labeled 
I A ( J ) . In Fig. 3, each numeral located rear a aot is the inoex I 
of that point. Each numeral locatec near the center of a line is 
the index J of that segment. 
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Fig. 3. Points and segments on a simple wire antenna 
lying in the xz plane. 


Fig. 4 shows the same antenna and the eight dipole modes 
defined by subroutine ISORT. The arrows indicate the reference 
directions for the mode currents and voltages. The mode index 
number I is placed near the terminal point 12(1). Mode I is a 
sinusoidal basis function which vanishes at the endpoints 11(1) 
and 13(1) and has unit current at the terminal point. These arc 
overlapping subsectional bases, and mode 1 extends over two inter- 
secting segments J A { I ) and JB(I). The reference direction for mode 
currents and voltages is from II to 12 toward 13. In Fig. 4, 
inodes 1 and 2 have terminals at the ground plane, with segment JA 
above and segment OB below the ground plane. This type of mode has 
no image. Modes 3, 4 and 5 have images. The size of the compressed 
matrix is NCM = 5. If we did not take advantage of the ground-plane 
symmetry, the matrix size would be N = 8. 


Table I presents some of the output data for this example 
and Table II lists the elements of the compressed matrix C(I,J) on 
return from subroutine IDANT. From Table I, the calculated imped- 
ance is Z 1 1 = 959 + j 664 ohms. For the same antenna with perfect 
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Fig. 4. Mode map for the antenna shown in Fig. 3. 


conducti vity, Zll = 879 + j 749 ohms. The calculated results 
should not be considered accurate without checking the behavior as 
the wires are subdivided into more segments. The longest segment 
should not exceed A/4. The thin -wire and delta-gap formulations 
are justified most readily if the wire radius does not exceed 0.007A 
Fortuitously, satisfactory results are often obtained for closed 
wire loops even when the wire radius is as large as 0.02 a. For 
dipoles, an upper limit of 0.007A is recommended. 

With 20-ohm resistive loads inserted in each end of each 
segment of the antenna in this example, the calculated impedance 
is 710 + j 206 ohms and the efficiency is 45.5 per cent with CMM = 1 
With and without lumped loads, the results obtained with the present 
computer program show satisfactory agreement with those obtained 
with the program in reference 1. (A new version of subroutine SCANT 
was used for these tests because the original version in reference 1 
does not always handle lumped loads properly.) 
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-.0305 A A 1 

.0031 183 

2 

.873 

91 . 1 

-.0000507 

.0027537 

3 

.739 

-77 . 8 

• O 3 3 A 892 

-.0022579 

A 

.271 

-3A .7 

• 03370AS 

-.0OOA880 

5 

.636 

-86.1 

.0001 362 

-.8320153 

EFF = 

93. AA Y 1 1 

= .003705 

- . 0 0i3A 88 

Zll* 959.07 66A.07 


PH 

TH 

D3P 

D’JT 

GPP 

GTT 

3 . 

85. 

O . 0 3 

.66 

3.03 

1.16 

20 . 

85 . 

-A 7 . A A 

1 • 2 A 

.00 

1 .33 

A 0 . 

85. 

- 3 9 . A 9 

2.6 3 

.00 

1 .83 

63. 

85. 

-33.35 

A . 1 7 

.00 

2.61 

80. 

85. 

-29. A6 

5.21 

.30 

3.A0 

100. 

85. 

-27 . 93 

5 .77 

.00 

3.77 
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r- co >o 


TABLE II 


Compressed Impedance Matrix 


I 

J 

C(I,J) 

C(0, 

I) 

1 

1 

16.1 

"j 

720.5 

16.1 

~j 

720.5 

1 

2 

7.6 

"j 

6.7 

7.6 

-j 

6.7 

1 

3 

-15.9 

"J 

1059.9 

- 7.9 

“j 

529.9 

1 

4 

-18.8 

-j 

16.8 

- 9.4 

“j 

8.4 

1 

5 

- 5.2 

+ j 

83.0 

- 2.6 

+ j 

41.5 

2 

2 

16.1 

"j 

720.5 

16.1 

“j 

720.5 

2 

3 

- 5.2 

+j 

83.0 

- 2.6 

+j 

41.5 

2 

4 

-12.7 

"J 

77.3 

- 6.4 

-j 

38.6 

2 

5 

-15.9 

"j 

1059.9 

- 7.9 

"j 

529.9 

3 

3 

21.2 

“j 

326.0 

21.2 

"j 

326.0 

3 

4 

- 9.2 

"j 

22.6 

- 9.2 

"j 

22.5 

3 

5 

- 9.2 

“j 

396.0 

- 9.2 

"j 

396.0 

4 

4 

51.3 

+J 

60.8 

51.3 

+j 

60.8 

4 

5 

22.1 

+ j 

420.9 

22.1 

+j 

420.9 

5 

5 

21.2 

-j 

326.0 

21.2 

"j 

326.0 


V. SUMMARY AND CONCLUSIONS 

A computer program is presented for a thin-wire antenna over 
a perfectly conducting ground plane of infinite extent. The analysis 
is performed in the frequency domain, and the exterior medium is free 
space. The antenna may have finite conductivity and lumped loads. 

The output data includes the current distribution, impedance, radia- 
tion efficiency and gain. The program uses sinusoidal bases and 
Galerkin's method and takes advantage of the ground-plane symmetry 
to reduce the storage requirements and computation costs. The 
subroutines are included in alphabetical order in the Appendices with 
a brief explanation. 
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APPENDIX 1. Subrout in? ANTI 


Subroutine ANTI is listed in Fig. 5. Between statements 14 
and 30, this subroutine sets up the excitation voltages CJ ( I ) and 
VJ(I) with the aid of the delta-gap model and the input data for the 
g«neratnr voltages VG(J). ANTI calls CR0UT to obtain a solution for 
the simultaneous linear equations . On return from CROUT, the dipole 
mode currents are stored in C0(1). The image currents are stored 
in CJ(K) in the DO LOOP ending with statement 80. The DO LOOP 
ending with statement 90 calculates the complex power input Yll and 
the time-average power input G. The power dissipated (DISS) in the 
lumped loads and the imperfectly conducting wire is obtained by 
calling PDISS. Finally, the radiation efficiency EFF is calculated. 

If IWP.CJ is positive, ANTI writes a list of the dipole mode 
currents CJ ( I ) . This list includes the normalized current magnitude, 
the phase in degrees, and the real and imaginary parts of the 
current. 


APPENDIX 2. Subroutine CROUT 


CROUT, listed in Fig. 6, solves a system of simultaneous linear 
equations with complex coefficients. This subroutine uses the method 
of P. D. Crout. Although this subroutine does net use pivoting, it 
is efficient and accurate in the present application. The input 
data are defined as follows: 


C(I,J) complex coefficients in the simultaneous equations 

S(I) excitation column 

ICC dimensions of C and S 

ISYM zero or one for symmetric or nonsymmetri c matrix 

IWR one or zero to write or suppress the solution 

112 one or two if C is original or auxiliary matrix 

N size of the square matrix C 


Of course, N must not exceed ICC. If IWR is a positive integer, the 
solution will be printed out with the following definitions: 


I index number of the solution S(l) 

SNOR normalized magnitude of S(I) 

SA absolute magnitude of S(l) 

PH phase of S(I) in degrees 

On the first call to CROUT, C ( I , J ) contains the original 
matrix. 112= 1 and CROUT generates the auxiliary square matrix, 
overlaying it in the same location C and destroying the original 
matrix, then CROUT proceeds to generate the solution, storing it in 
S ( I ) and destroying the original excitation coiurin. 
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Next we might wan* another solution of the same system of 
simultaneous linear equations hut with a new excitation column. This 
coulji be obtained by recalculating the original matrix C(I,J) and the 
new excitation column and calling GROUT again with 112 = 1. However, 
there is no need to recalculate C(I,J). Instead generate the new 
excitation colunn, set 112 = 2 (or any integer other than 1) and call 
CROUT again. GROUT uses less cocput.er time when 112 differs from 1. 


APPENDIX 3. Subroutine EXPJ 

Subroutine EXPJ, listed in Fig. 7, evaluates the exponential 
integral defined as follows: 


V 2 


W12 = 


— -- ^ v = E 3 (VI ) - E 3 (V2) + j 2mr 
v 


VI 


where the integration path is the straight line from VI to V 2 on the 
complex v plane and 


ct> 



z 


Ej(z) denotes the principal branch of the exponential integral. To 
generate W12, EXPJ calculates E 3 { V 1 ) , subtracts E 3 ( V2 ) and adds j2nTr. 
The integer n is zero unless the straight-line integration path 
intersects the negative real v axis at a point between VI and V2. 

When there is such an into, section, n s 1 if VI lies in quadrant 1 or 
2 and n = -1 if VI lies in quadrant 3 or 4. 


APPENDIX 4. Subroutine IDANT 


Subroutine IDANT is listed in Fig. 8. This subroutine stores 
the quantities CDK(J) = cos kds and SDK(J) = sin kdj where dj is the 
length of seqment j. The proqram wri tes AK, DMAX arid DM I N and aborts 
if 


a. the length of the shortest segment is less than the wire 
radius, or 

b. the longest segment has a length d such that kd exceeds 3, 
or 

c. the wire radius a is such that ka exceeds 0.1. 

IDANT calculates the elements in the compressed inpedance 
matrix C(I,J) as follows. Select a source segment K and a receiving 
segment L, where K and L range from 1 to NS. The mutual impedances 
PI 1 » P12, P21 and P22 between the two segments are obtained by calling 
ZGMM if K = L, ZGMM if the segments intersect, or ZGS if segments K 
sr.o L do not intersect. 

Now select a test dipole I sharing segment K, and an expansion 
dipole J sharing segment L. Add the appropriate segment-to-seginent 
impedance to the dipole-to-dipole impedance C(I,J). When this 
procedure has been completed at statement zOO, the impedances C(1,J) 
are appropriate for a perfectly conducting thin-wire system with no 
lumped loads. 

Between statements 200 and 262, the impedance matrix C(I,J) is 
modified to account for the finite conductivity of the wire antenna. 

The surface impedance ZS is obtained by calling ZS'JRr. For each 
segment K, the program selects a test dipole I and an expansion dipole 
J sharing this segment. The contribution to C ( I , J ) associated with 
finite conductivity is ZSAM if dipoles I and J have terminals at the 
same end of seqment K, and ZOPP if they have terminals at ooposite 
ends. C(I,J) is not affected unless dipoles 1 and J share one or 
two segments. 

Between statements 262 and 280, the i npe dance matrix is nodified 
to account for the lumped loads. Each diagonal element C(I,I) is 
modified by adding the impedance of the lumped load inserted at the 
terminals of rode I. If modes I and J share a segment and have termi- 
nals at the same point, C ( I , J ) is modified by adding or subtracting 
the impedance of the lumped load inserted at the terminal end of 
this segment. (Add or subtract ZLD if node currents 1 and J have the 
same or opposite reference directions on the shared segment.) 
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APPENDIX 5. Subroutine IFFLD 


Subroutine IFFLD, listed ir. Fig. 9, calculates the far-zone 
field of the thin-wire antenna. 

Let (r,0 ,<) denote the spherical coordinates of the distant 
observer, and let E 0 (I) and E.(I) denote the elc-ctric field intensi- 
ties of dipole mode I with unit terminal current. Then 

EPP(I) = (r/x) eJ’ kr E t (l) 

En(I) = (rA) c-* kr E e (I) 

The field of sinusoidal dipole mode I may he regarded as the sum of 
the fields of each of its two segments. The field of segment K is 
obtained by calling subroutine ZFF, and EPP(I) and ETT(I) are 
generated by adding the appropriate numbers obtained from two 
different calls to ZFF. In the [X) LOOP ending with statement 260, 
the antenna field is calculated as a weighted sum of the mode fields 


CJ(I) EPP(I) 


CJ(I) ETT(I) 


where CJ(I) denotes the terminal current of rode I and EPH and ETH 
denote the dimensionless range -in depen dent form of the antenna fields 
E; and E 0 . 

G denotes the time-average input power to the antenna, and 
GPP and GTT are the < -polarized and O-polarized power gains. Sub- 
routine IFFLD is called once for each angular direction. In the 
input data supplied to this subroutine, PH and TH denote t a"d 0 in 
degrees. 

This subroutine is useful for wire antennas with or without a 
ground plane. For an antenna over a ground plane, IFFLD must be 
supplied with information on the complete system including the image. 


as follows: 


N 

EPH = 1 

1 


N 


ETH = 


16 




APPENDIX 6. Subroutine ISORT 

Subroutine ISORT, listed in Fig, 10, is described briefly in 
Section III. This subroutine sets up the image segments and points 
and calculates the segment lengths. Then it checks the input data 
for consistence. The data are considered inconsistent and the run 
is aborted i f 

a-. NPGP is greater than zero but one of the points on the 
ground plane has no segment (with index J less than or 
equal to NSGP) connected to it, or 

b. a real point situated above the ground plane has no 
segment (wi tli index J less than or equal to NRS) 
connected to it. 

Between statements 32 and 50, this program calculates the 
number of modes N on the complete structure. The run is aborted if 
the dimensions are inadequate. 

Between statements 50 and 58, the program sets up the modes 
that will not have images. The number of modes of this type is 
NPGP, and these modes have the lowest index numbers. Hode I has 
terminal point 12(1) = I on the ground plane, endpoint 11(1) is above 
the ground plane, endpoint 13(1) is the corresponding image point 
below the ground plane, and segment JA(I) is the lowest -numbered 
real segment with endpoint I. 

Between statements 58 and 65, the program sets up the rest of 
the real modes, (lodes of this type have the terminal point 12 on 
or above the ground plane. Each of these real modes (with index I 
greater than NPGP) has an image which is established between 
statements 65 and 75. 

Below statement 75, the last part of the program counts the 
number of dipole nodes sharing segment J, denoted by N0(J). It also 
stores a list of the dipole modes sharing segment J, denoted by 
MD(J,K). A segment may be shared by as many as four modes. 


APPENDIX 7. Subroutine PDISS 

Subroutine PDISS is listed in Fig. 11. This subroutine 
calculates the time-average power (DISS) dissipated in the lumped 
loads and the imperfectly conducting wire. The power is calculated 
for one segment at a tine, and the total power dissipated is the sum 
of the powers dissipated on the various segments. On segment K, CJA 
and CJB denote the currents at endpoints I A ( K ) and IB(K). RLA and 
RLB denote the lumped resistors inserted in segment K at endpoints 
IA and IB. 
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This subroutine is suitable for a wire structure in free space 
and also for a wire structure over a perfect ground plane. If there 
is no ground plane the total nunber of segments NS must be supplied 
as the tenth calling parameter, and D1SS denotes the power dissipated 
on the entire structure. If there is a ground plane, the number of 
real segr.. ,its MRS is supplied instead, and DiSS denotes the power 
dissipated on the rea. segments. 


APPENDIX 8. Subroutine ZFF 

Subroutine ZFF, listed in Fig. 12, calculates *he far-zone 
field of a sinusoidal electric nonopole in free space. The monopole 
has endpoints at (XA.YA.ZA) and (X'5,YB,ZB).. (These symbols denote 
kx, ky and kz.) Let E denote the electric f i e 1 d intensity. The 
dimensionless range -in depen dent *ield is defined by 

F * (r/>.) e jkr E 

EP1 and ET) denote F, and F 0 for the moae with unit current at 
(XA,YA,ZA). EP2 and’ET2 denote F 4 and F r , for the mode with unit 
current at (XB,YB,Z8). The far field vanishes in the endfire direction 
where GK - 0. 


APPENDIX 9. Subroutine ZGMM 

Subroutine ZGMM, listed in Fiq. 13, calculates the mutual 
impedance between two f i 1 amenta ry monopoles with sinusoidal current 
distributions. The dipole-dipole mutual impedance in Eq. 20 of 
reference 2 is the sum of four monopole -monopole mutual impedances. 

The monopole impedances are calculated by ZGS with Simpson's rule or 
by ZGMM with closed-form expressions in terms of exponential integrals. 

If the monopoles are parallel, let the z axis be parallel with 
both monopoles. The coordinate origin may be selected arbitrarily. 

SI and S2 denote the z oordinates of the endpoints of the test 
monopole, T1 and T2 are the z coordinates of the endpoints of the 
expansion monopole, and D is the perpendicular distance (displacement) 
between the monopoles. The mutual impedance of parallel monopoles is 
calculated in the last part of ZGMM below statement 110. 

For skew monopoles, let the test ronopoie s lie in the xy 
plane and the expansion monopole t in the plane z = D. (D is the 
perpendicular distance between the parallel planes.) If the monopoles 
are viewed along a line of sight parallel with the z *xis, the ex- 
tended axes of the two monopoles will appear to intersect at a point 
on the xy plane. Let s measure the distance along the axis of the 
test monopole with origin at the apparent intersection. SI and S2 
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denote the $ coordinates of the endpoints of the test monopole. 
Similarly, let t measure distance along the axis of the expansion 
monopole with origin at trie apparent intersection. T1 and T2 denote 
the t coordinates of the endpoints of the expansion moncpole. Let 
s and t be unit vectors parallej with the positive s and t axes, 
respectively. Then CPSI = s • t = cos «•*. The monopole lengths are 
d s and dt, and the remaining input data are defined as foi'ows: 

CGDS = cos kd s 

SGD1 = sin kd s 

SGD2 - sin kdt 

ZGMM calls FXPJ for the exponential integrals. ZGMM is 
specialized for sinusoidal monopoles in free space. In ZGMM the 
input data SI, S2, Tl, T2 and D denote ks ]f ks^» kt 2 , kt-> and kd, 
respectively. Otherwise, ZGMM is the same as GGMM.‘ 

The output data from ZGMM are the impedances PH, P12, P21 
and P22. In defining these impedances, the reference direction is 
from SI to S2 for the current on monopole s, and from Tl to T2 for 
the current on monopole t. In the impedance Pij, the first subscript 
is 1 or 2 if the test dipole has terminals at Si or S2 on rnonopole s. 
The second subscript is 1 or 2 if the e*pansion dirx)le has terminals 
at Tl or T2 on monopole t. The endpoint coordinates SI, S2, Tl and 
T2 may be positive or negative. The monopole lengths d s and d^ d'e 
assured positive in defining the input data CGDS, SGD1 and SGD2. 

For parallel monopoles, CPSI = 1 or -1. SI, S2, Tl and T? 
are cartesian coordinates for parallel monopoles and spherical 
coordinates for skew monopoles. For skew monopoles, the radial 
coordinates SI, S2, Tl and T2 tend to infinity as the angle •; tends 
to zero or r . Therefore, if the monopoles ire within 4.6 degrees 
of being parallel, they are approximated by parallel dipoles. 


APPENDIX 10. Subroutine ZGS 

Subroutine ZGS, listed in Fig. 14, calculates the mutual 
impedance between two filamentary monopoles with sinusoidal current 
distributions. (The dipole-dipole mutual impedance in Eq. 20 of 
reference 2 is the sum of four monopole -monopole mutual impedances.) 
The endpoints of the axial test monopole s an? at (XA,YA,ZA) and 
(XB,YB,ZB), and the endpoints of the expansion monopole t a re at 
( X 1 , Y I , Z 1 ) and (X2,Y2,Z2). DS and DT denote the lengths of nono- 
poles s and t. Dimensionless forms are used for the input data. 

For example, XA, AK, DS and DT denote kx^, ka, kd$ and kd^. CAS, 

CBS and CGS are the direction cosines of monopole s, and CA, CB and 
CG are the direction cosines of nonopole t. 
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If INT = 0, 2GS calls ZGMM for the closed-form impedance 
calculations. Otherwise ZGS calculates the mutual impedance via 
Simpson‘s-rule integration with the following number of sample 
points: IP - INT + 1." If the rnonopoles are parallel with small 

displacement, ZGS calls ZGMM to avoid the difficulties of numerical 
integration . 

For the fields of the test nonooole, ZGS uses Eqs. 75 and 76 
of reference 2. The current distribution on the expansion nonopole 
is given by Eq. 74 of reference 2 . With an origin at (X1,Y1,Z1), 
the coordinate T measures distance along the expansion monopole. 

Thus T is the integration variable. 

Let the coordinate s measure distance along the test monopole 
with origin at (XA,YA,ZA). From any point T on monopole t, construct 
a line to the test monopole such that the line is perpendicular to 
the test monopole. SZ denotes the s coordinate of the intersection 
of this line with the test monopole. The length of the line is the 
radial coordinate o, and RS denotes r 2 . R1 and R 2 are the distances 
from (XA,YA,ZA) and (XB,Yb*ZB) to the point T. Cl is the current 
at T for the mode with terminals at (X1,Y1,Z1), and C2 is Ine current 
at T for the other node with terminals at (X 2,Y2,Z2). C denotes 
the Simpson 's-rule weighting coefficient. 

Below statement 300, ZGS performs some analytic geometry in 
preparation for calling ZGMM. The remaining part of this Appendix 
concerns this last part of subroutine ZGS. 

Let s denote a unit vector from (XA,YA,ZA) toward (XB,YB,ZB), 
and let ? A denote a unit vector from (X1,Y1,Z1) ta/ard (X2,Y2,72). 

Tnen s • t = cos v = CC where v is the angle formed by the axes /f 
tne two rnonopoles. Let monopole s lie in one plane P s and rnjnopole t 
in another parallel plane P^. £AD, CBD and CGD are the direction 
cosines of the unit vector d = t x s / sin v which is perpendicular 
to botn planes. To obtain the distance DK between the planes, we 
construct a .vector R. . from (XA,YA,ZA) to (X1,Y1,Z1) and take 
DK = R n • d. 

Construct a line from (X1,Y1,Z1) to the test monopole, such 
that the line is perpondicul ar to the test monopole. SZ denotes the 
s coordinate of the intersection of this line with the test monopole, 
and the cartesian coordinates of this intersection are XZ, YZ and 11 . 
The direction cosines of $ x d are CAP, CBP and CGP. 

From the point (X1,Y1,Z1) in plane construct a perpendicu- 
lar line to the point (XP1 f YP1 ,ZP1 ) in the plane P s . This line is 
parallel with cj and has length DK. Let R represent a vector from 
fXZ,YZ,ZZ) to (XP1,YP1,ZP1). PI denotes R-(s x a). SI and T1 are 
defined in Appendix 9. 
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Subroutine ZGS is essentially the same as GGS except the 
medium is specialized to free space in ZGS. 


APPENDIX 11. Subroutine ZSURF 

Subroutine ZSURF, listed in Fiq. 15, calculates the surface 
impedance of a solid circular-cylindrical wire with exterior exci- 
tation. ZS denotes the surface impedance in ohms, and the input 
data are defined as follows: 

AK ka, where a is the wire radius 

CMM conductivity of the wire in megamhos/m 

FMC frequency, MHz 

The surface impedance is defined by Z s - E^/H. when* the fields E and 
H are evaluated at the surface of the^wire. Yhis subroutine calcu- 
lates the impedance *or the lowest order cylindrical mode with fields 
Z 7 and H. independent of <. The wire is considered to be a qood 
conductor in the sense that the displacement current is negligible 
in comparison with the conduction current. In the present applica- 
tion, we require the surface impedance appropriate for the current 
distribution I(z) * sin kz. For a highly conducting wire, however, 
this impedance is considered to be the same as that for a urifomly 
distributed current. BER, BE I » BERP and BE I P denote the Kelvin 
functions her, bei, ber' and bei' with arqument x. When x is less 
than 8, BES and BES1 denote the Bessel functions J 0 and J, w.th 
argument 

z = x e~ jTT/ ^ • 


When x is greater than 8, 


SUHftnuT INF ANT J < J A, 1B,-I1 , 12. I 3, iWKCJ, ltfR J1 F . JJ 2 f ICC, JNS, JA, JH, 001 

?JPP,MO.N,NCM,Ni;,N|.h,NPGP, Nk S » NS • C , CPK « Sl»K , C J « Cmm , 0 . k F F , G r VG • V.) , 002 

3Y1 1, 7 11 1 / H, 2 Li) ) 003 

CUNPLF x C( ICC. Iff > .CJC 1 ) * Vf,( 1 ) , VJ( 1 ) , ZLDI 1) 004 

CMMPLFX Yll,211,/M,r.JI, VJI.PYY 006 

DIMENSION 1 A ( 1 ) , I *i ( 1 ) , iO ( 1 ) ,COK( 1 ) , SOX ( 1 ) • I'M 1 i . I 1 ( 1 ) * 1 2 ( 1 ) * I 3 ( 1 ) 006 

D1HFNSICW Ml)( INS^K JA( 1 ) • Jh( 1 ) , I Ok KM 1 ) , JGLM 1 : 007 

1 FORMATCHX, Mr « f t)X, ‘Vb( J > = ' , 2F 1 <». 2 ) 000 

2 FORMAT I i.x , »js • . | ■> , 6X, • 7 l Of J ) r • f ?» 1 0. ? ) 009 

3 FORMAT ( 1 OX , « I * , AX , ' :»AGm! I UOH HX, • PMASb ' » 9X , • kE AL ' , 8X, • 1 '-iAG I N/.R Y 1 > 010 

4 FORHAUh.HlO.lfMO.BilFlO. l,?Fb.7| Oil 

6 K1RMA1(1h.)» OIF 

1 F ( I 1 2. 0 r . 0 ) 112-1 013 

I F ( IWR11F.LF.01 GO 1 U 14 01 A 

nn io j= i , nrs 016 

A VG - C. APS ( VC* ( J ) ) 016 

10 IF( AVG.GT.O.Ol )WRI1E(6, 1 )J«VG(J) 017 

W R 1 T F ( 6 » 6 ) OIK 

iPiNin.u . oioo in !'♦ 019 

7MAX-.0 020 

on 12 J= 1 * MR S 021 

A 2 L =C. A H S T 7 L 0 ( J) ) 02? 

1FIA7L.LT.0. 1)01) T;3 12 02 3 

2 M A X = 1 • 024 

W R 1 T F ( 6 , 2 ) .) , / L I M J ) 026 

1? CONTINUE 026 

IF ( 7MAX.GT. 0. 6 )WR I 1 k ( 6, 6 ) 027 

14 no 30 l = 1 , N C M 02 8 

CJI I ) = (.O f .0) 029 

L =2 030 

F«C~ 1 • 0 :> i 

K = J A ( 1 ) 032 

IF( 1.GT.NPGP1GU TO 15 033 

L = 1 034 

F AC = 2. 035 

I F ( J rt ( 1 ) .LT.K ) k = Jrt ( I ) 036 

15 00 26 *K=1,L 037 

K A = 1 A ( K ) 038 

K B = I B ( K ) 039 

J J = K 040 

F 1 = F AC 041 

1 F ( KB. EO. I 2 ( I ) ) GO TO 22 04? 

JF( KH. EO. I) { l))h I --FAC 043 

cjiihCjm+FMvcuji 044 

GO TO 26 045 

2? IF <KA. FO. I 31 I ) )F !=-FAC 046 

JJ=K4KRS 047 

CJ( 1 )=CJ( 1 >*F1«VG( JJ) 048 

25 K = .) B ( I ) 049 

VJ(IXJII) 060 

K=I+JPP 051 

30 I F ( l • G T ■ N P GP ) V J (K ) = - V J ( I ) 062 

I SYM= 1 053 

1FCN.F0.NPGP) ISYM=0 064 

I W R = 0 065 

CALL CROUT(C*CJ, ICC, 1 SYM f 1WR, I 12»NCM) 066 

112=1? 067 

CMAX-.O 058 

on 80 U1,NCM 059 

C A =C AH S ( C J ( 1 ) ) 060 

K=I 4 JPP 061 

IF( I.GT.NPGP)CJ(K) =-CJ( I ) 062 


Fig. 5 . 
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Subroutine ANTI. 

22 


80 lb ( CA.GT. CMAX>C«AX = CA 063 

l»-( IWRCJ.GL. 1 )W*ITI (6, 3) 064 

Yll-(.0».0) 066 

G=1 . 066 

ff-F=.0 067 

niM.Oi.O) * 068 

IF (CMAX. LF . 0. >00 TO 600 069 

Ca.O 070 

bn 90 I = 1 , N 071 

CJlsf.JUl 072 

VJ I =VJ( n 073 

OYYaC J 1*C0NJG( V.ll 1 074 

lb I I . I f-.OO!) VI r- Y1 1 ♦ DY Y 076 

G = C»4KF Al. ( OYY ) 076 

IF (I wRf.J. IF. 0)(,n 111 90 077 

If ( 1.01. NCM)Gll TO 90 078 

CA = CA8S(CJ1 ) /CMAX .. 079 

PH-.O 080 

IMCA.CT. I . f - 30) PH = 6 7. 2967MV AT AN?( A I MAGI CJ 1 ) «KbAL (C J t ) ) 081 

Wk I T F 1 6, 4 ) 1 ♦(. A, PH, C J 1 082 

90 CONTINUE 083 

IN IwkCJ.GE. 1 )UK1T| (6, 6) 084 

G=G/2. 08 8 

711=1. /Yll 086 

f F f s 1 00. 08 7 

1FTCMM.LF.0. . A Ml). NLD.LE.OIGO TO 800 0R8 

CALL Pill SS ( 1 A, I B t INS» 11. I 2 » I 3* Ml), KD, ML D» NK 5*CJ*CMM,D» COK » 089 

?snK«niss» 7 H v uiM 090 

F H r 1 00. * ( G- 0 1 SS 1 /G 09 1 

800. RF TURN 092 

FNO 093 


Fig. 5. Subroutine ANTI - continued 



SUBROUT I NF CROU T ( C » S » ICC, iSYM, J WR f j \2 r N) 

001 


COMPLEX f ( ICC, ICC). S( n 

002 


COMPLEX F,P,SS.l 

003 

2 

FORMAT! IX, 1 IS, IF 10. 3, 1 FIS, 7, IF 10. 0) 

OOA 

5 

FORMAT! !►«/)) 

005 


1 F ( 1 12. NF . 1 ) GO TO 22 

006 


I F ( N , F 0 . 1 )S( 1 ) ~ S ( 1 )/C( 1,1) 

007 


1F( N.t 0.1)GO 10 100 

008 


IF! 1 SYM. 98.0)00 10 h 

009 


00 8 I = 1 » N 

010 


DO 8 J=>l , N 

Oil 

6 

C ! J, I 1 = C( 1 , J 1 

01? 

R 

F=C<1 , 1 ) 

013 


DO 10 l=?,N 

01 A 

10 

C< 1,1 >=C!1 ,L ) /F 

CIS 


DO 20 1=2, N 

018 


LLl =L-1 

01 7 


00 ?0 I r L » N 

018 


f - c ( i , ) 

019 


DO 11 K= 1 , L L L 

020 

1 1 

F = F - C ! 1,K)»C(K,L) 

021. 


C ( 1 , L ) = F 

02? 


IFIL.EO.DGO 10 20 

02 3 


P=C(l ,u 

02 A 


IF! lSYM.EO.OJGO 10 IS 

025 


F «C ! L, 1 ) 

026 


00 1? K=1.LIL 

021 

1? 

FrF-C.I L , K ) * C ( K , 1 ) 

028 


C(l , 1 )=F/P 

029 


GO TO 20 

030 

1 5 

F=C( I ,L ) 

031 


C ( L , I ) =F / P 

03? 

20 

CONI INOE 

033 

22 

DO 30 l = 1 , N 

03A 


P = C i L , L ) 

035 


T = S ( L ) 

036 


1 F ( L . FO. 1 )G0 TO 30 

037 


UL=L-1 

03R 


DO 2S K= 1 , l L L 

039 

2b 

T=T-C(ltK )*S(K ) 

0 AO 

30 

S ( L i = T / P 

OA 1 


DO 38 1=2, N 

0A2 


1 =N-l -*• 1 

OA 3 


11=1+1 

0 A A 


T = S( 1 ) 

0A5 


o 

uf 

* 

II 

OAO 

35 

7 rT-C 1 l t K)* S!K) 

OA 7 

38 

S( 1 ) = T 

0 A 8 


IF ( lwR.tr . 0 ) GO TO 100 

OA 9 


WR 1 1 EC 8, S ) 

050 


C NOR = • 0 

0S1 


DC) AO 1 = l , M 

05? 


SA=C ABS ( S ( I > ) 

093 

'♦0 

IF ( SA. GT.CNDR )CVOR=SA 

05A 


1 F ( CNOR. LF .0. )CN0R~1. 

OSS 


DO A A 1 = 1 , N 

056 


SS = S( I ) 

05 7 


$ A = C A {) S ( S S ) 

058 


SNOR=S A/CNOR 

059 


PH = . 0 

060 


If ! SA.G1 .0. )PH = S7. 295 7 8* A TAN? ( A I MAG ( SS ) • RE AL ! SS ) ) 

061 

AA 

WR1 1 E ! 8, 2 ) I , SNlM , SA , PH 

062 


WRITE (8, S) 

06? 

1 00 

RETURN 

06A 


END 

065 


Fig. 6. Subroutine CROITT 
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SUBROUTINE EXPJ( VUV2, W12) 001 

COMPLEX EC«E1*> V S.T.UCf VC* V1»V2* W12.2 002 

DIMENSION V( 21 ), W( 21 ),U( 16),6( 16) 003 

DATA V/ 0.222846676 OG , 004 

20. 118893216 01,0.299273636 0), 0.577514366 01 , 0. 983746748 01, 009 

20.159828748 02 , 0.9 33078 1 26-01 f 0. 49269 1 74b 00 , 0 . 1 2 1 559548 01. 006 

20. 226994958 01,0.366762271: 01,0.542933666 01,0. 796991628 01, 007 

20.1012022HE 02, 0. 1 3130232F 02,0.166966088 02,0.207764796 02, 008 

20. 286238948 02, 0. 3] '.079 1 9c 02 , 0. 3893068 3E 02 , 0. 4H0260B6E 02/ 009 

DATA W/ 0.468964608 00, 010 

20. 1 7000831- 00,0. 1 13373388 00, 0. 1 0399 1 978-01 ,0. 261 01 7206-03, Oil 

20.89894791 6-06, 0. 2 18236878 00,0.342210176 00 , 0. 263027986 00, 012 

2 0. 126429 826 00, 0. 402 06865E-0 1 , 0. 89*38 7 l HE -02 , 0. 1 2 1 2'. 36 1 6-02, 013 

20. 1 1 1 6 7440b -03 , 0. 6 4 999267 E -09 * 0. 22263 1 69 E -06 , 0 • **22 7*3046-08, 014 

20. 392189738-10,0. 149691928-12,0. 148302/06-2 9, 0. 160099496-19/ 015 

DATA D/ 0.224WSR426 02, 016 

2 0.744119688 02,-0.414319768 03,-0.787543398 02, 0.112547448 02, 017 

2 0.160217618 03,-0.238621996 03 , -0. 500<>46H 7 E 03 • -0 . 68487 8548 02, 018 

2 0. 122547786 02, -0. 1 0 1 6 1 9768 02, -0. 472 1 959 1 6 01, 0. 79729681c 01, 019 

2-0.210695748 02, 0.220464906 01, 0.89728244b* 01/ 020 

0 A 1 A 6/ 0. 21 1031076 02, 021 

2-0. 379597ft H 03,-0.974892208 02, 0.129006728 03, 0.179492266 02, 022 

2-0.129109316 03,-0.957095746 03, 0.136248018 02, 0.146967218 03, 023 

2 0.179495288 02,-0. 329810148 00, 0.310288366 02, 0.H1657697E 01, 024 

2 0.222369618 02, 0.391248926 02t 0.816367996 01/ 029 

2 = V 1 026 

Of) 100 J I Ms 1,2 027 

X = R 8 A 8 ( Z ) 028 

Y = A l MAG ( Z I 029 

615-1.0*. 01 030 

AB = C AOS ( Z ) 031 

1 8 ( AR. 60.0. )f,n T(1 90 032 

)F(X.G6.0. .AND. AB.GT.10. ) GO TCI 80 033 

YA = ABS ( Y ) 034 

1 F ( X . L 6 . 0 . .AND. YA. GT. 10. IGO ID 80 035 

1 F ( YA-X.GE . 1 7.5.0R. YA.GE. 6. 5. UR. X4YA.GE. 5. 5.0R.X.G6 .3. 1G0 TO 20 036 

JFlX.ie.-'?. IGO TO 40 037 

IF( YA-X.G6.2. 5) GO 1U 90 038 

!6< X4YA.G6. 1.9IG0 10 30 039 

10 N = 6. -*-3 . * AB 040 

615=1. /(N-l. )-Z/N**2 041 

15 N = N — 1 042 

FI 5=1. /(N-l. 1-7*815 /N 043 

1F1N.G6.31GO TO 15 044 

8 15 = Z*6 15-CMf'LXl . 5/7216+AL OG ( A B ) , AT AN2( Y, X ) 1 045 

GO 1(1 90 046 

20 Jl=l 047 

J2 = 6 048 

GO TO 31 049 

30 J 1 = 7 060 

J2”2 i 051 

31 S=(.0,.0) 092 

Y S = Y * Y 053 

DO 32 I = J 1 , J2 054 

X 1 =V( I) *X 055 

CF=W( 1 ) / ( X I * X 1 ♦ YS ) 056 

32 S = S^CMPLX ( X 1 *CF, -YA*CF ) 057 

GO TO 54 OSH 

40 T3=X*X-Y*Y 059 

T4=2.*XCYA 060 

T5=X*T3-YA* 14 061 

1 6=X*1 4+YA* 1 3 062 

Fig. 7. Subroutine EXPJ. 
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UC=CMPLX(f)( 1 1 )*0< 1? )*X*r»t i:4)*T3*T W< 12 )*YA-F( I *3 ) «= T a . 06 3 

2 E( 11 ) + F( 1?>*X+E( 13)*T3 + 1(>*0(12)*YA + {J< !3)*14) 06A 

VCrCMPLX(D( 14)4U( 15jcX+g( 16RT34T5-F( 15 R V4-c( 16 1* I A, 065 

2 E ( R ) + fcf lb )*X*£| lb)«T3*Tb+l)( l*>)*YA4f)(16) « 14) 066 

GfJ 1(1 52 067 

50 T 3 = X* X-Yo Y 068 

TA=2.*X*YA 069 

T5=X*T 3-YA*T4 070 

T6=X* TA+YA* 7 3 . 071 

T7 = X*T6 ■ Y A $ T 6 072 

1 8 = X*> 1 6 + Y A* (5 073 

T9=X*T 7-YA*l 8 07 A 

T 1 OsX* 1 «-» YA* T7 075 

UC=CMPLX(D( 1 ) 4i)f ? )*X4i)( 3 )«1 3*U( 4»*T54|).| 6RT 74 T9-( E( 7RYA4E ( 3) *T4 076 

24ER)*T64FC6RT8) f MlRE(2)*X4E(3l« : T34E(4)*Tb4EC6Rl74T104 077 

3<0(2l*YA4l>( 3)*R4l)RRT6 + |)| 5)* T8)) 07 6 

Vt=CMPL X( l)( 6 )♦()( /)».U0(B)#T3*0(<1)CT^* >! 1 0 )<* T 7 4 T Vr ( F ( 7 ) * Y A ♦ E ( 8 ) * 7 A 079 

2 + E (9 RT64E( 10 1* lfi) ,F(6)4E( /RX4E(6)*T34t(9RT54E( 10) *1747104 080 

3(0( 7 RYA4f)( 8 )*TA4D( 9 )*T64|)( 10 RTft 1 I 081 

52 FOUC/VC 082 

S-l C/CMPLX( X, Y A ) 083 

5A FX=FXP(-X) 08 A 

T=Fx*CMPt_X(C05>( YA I, -SI? V (YAM 08 5 

E 1 5 = S * T 086 

56 IF ( Y. t T.O. )fc lb*CONjG( E15 ) 087 

GO 10 90 088 

8 0 E 15c. 4093 19/ ( 7 ♦ . 19 30 A A )4.4?]&31/{24]. 02666 ) ♦ . 1 A 7 1 26/ ( 7 4 2. 56 788 ) ♦ 089 

2.206 335E-1 /( Z«A .90C35 M. 1 0 7A« i t -2/ ( 7 *8 . 162 15)4. 1 5865Afc-A/ ( /♦ 090 

312. 7342) '•3170316-7/(2419.3957) 091 

F 1 S=F 1 5*CEXP( -/ ) 092 

90 ir ( JIH.cO. 1 1W12-616 093 

100 7 =V? 09 A 

Z S V?/V l 095 

1 H = A 1 A N 2 ( A 1 M A G ( l ) , PEAL (Z ) ) - A T AN2 ( A 1 H AG ( V2 1 *KEAL ( V2) ) 096 

24 AT AN?( A I MA(>( V 1 ) , REAL ( VI ) ) 09 7 

A tt = A H S ( 1 H ) 098 

If (AH. LI. 1. ITHr. 0 099 

I F < TH.GT. 1 . 1 TH=6. 2631853 100 

I F ( TH. L T. - 1 . 1 lH=-6. 2831 853 101 

W12-W12-E154CMPLX( • 0 » 1 H ) 102 

RETURN 103 

END 1 0 A 


Fig. 7. Subroutine EXPJ - continued 
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SUBROUT INF lOANTt ! A f IB, ICC, INS. I NT, I ] . 1?, I 3 , J A , JB , JPP , MU, N, NCM, 00 1 

2ND, NLD,NP, NPGP, NRS ,NS, AK t C, CMM, D, FmC,CI)K, SDK f X , Y , 7 , / H, 7ll)) 002 

COMPLEX 7. SfZH f P(2.2)»0(2*2)*ClJtiSAh ? Z0PF 003 

COMPLEX C( ICC* ICC ) «2L1)( i ) 00* 

DIMENSION X( 1 ) f Y( I) , 7 ( 1 ) . I A( l ) , IBM ) « N1H ) ) , CDK ( 1 1 , SDK ( 1 ) , D( 1 ) 005 

DIMENSION I 1 1 1 ) , I 2 11) , I3t 1 )• JAt 1) , JBt 1 ) , MUt INS, *1 006 

DATA TP/6. 28 3 1H/ 007 

2 FORMAT ( 8X t • AK = • f F H .6, 5X f • UMAX* ' , F8.* , 5X, H)MlN= • , F8.* ) 008 

DO 10 1 = 1, NC M 009 

DO 10 J = 1 , NC M 010 

10 C( I , J ) s ( • 0, • 0) Oil 

DMAXa.o 012 

DMJN=200. 01 3 

DO 20 J * 1 , NR S 01* 

DJ *D( J ) 015 

IF (DJ.GT.DMAX JUMAX^DJ 016 

IF(DJ.LT.DM1MDM|N--DJ 017 

CDK ( J ) =COS I DJ ) 018 

SDK ( J ) * S I N ( D J ) 019 

K*J*NR$ 020 

CDKI K ) =CDK( J ) 021 

20 SDK ( K ) s SDK ( J ) 022 

IFtDMIN.LT. AK)G(! TO 21 023 

1 F < DMAX.GT . 3. )f.n TO 21 02* . 

I F ( AK.GT.O. I ) GO TO 21 025 

GO TO 22 026 

21 WR ! T6 ( 6» 2 ) AK, DHAX » DM I N 027 

N=*0 028 

RETURN 029 

22 DO 2 00 K « 1 , NS 030 

NDK-NDi K ) 031 

KA«]A(K) 032 

KB-IB(K) 033 

DK*D(K) 03* 

DO 200 L « 1 , NS 035 

NDL«=ND( L 1 036 

LA* 1 A( l ) 037 

l B * 1 B t L ) 038 

DL «=D ( L ) 039 

N I L c 0 0*0 

DO 200 II -1, NDK C*1 

I *MD t K , I I ) 0*2 

I F ( 1 . G T • N C M ) GO TO 200 0*3 

F 1 =1 . 0** 

IFtKB.FQ. I2t I ) >GG T Cl 36 0*5 

IFtKB.tO. 1 1 t l) ) F l * - 1 • 0*6 

I S = 1 0*7 

GO TO *0 0*8 

36 IMKAJU.I3II))M=“1. 0*9 

1 S c 2 050 

*0 DO 200 J J = 1 , NDL 051 

JsMDIlfJJ) 052 

I F ( I • GT* J )G0 TO 200 053 

F J * l . 05* 

IFtLB.FO. I2t J) )G0 TO *6 055 

IF(LH.tU.IlIJ) ) F J * - 1 . 056 

JS* 1 057 

GO TO 50 058 

*6 IFtLA.FO. I3(J) ) F J * - 1 • 059 

J S *2 060 

50 IFtNlL.NF.OlGO TO 168 061 

N I L = 1 062 


Fig. 8. Subroutine IDANT. 
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IFfK.bO.L )Gll TO 1 20 063 

1ND={ L A-KA )* < LB-KA )*( LA-KB )* ( LH-KH I 064 

IF ( INO. FO.OIGO l’u 60 066 

C SFGMEN1 S K ANO L SHAKE NO POINTS 066 

CALL 7C. S(X(KA)*Y(KA) # Z(KA),X(KH)*Y(Kft)*Z(KB)*X(LA)*Y(LA)*Z(LA)f 06 7 

?X( LB ) , Y( LB) , / (LH) , AK « t)K * C OK ( K ) , SDK( K ) , DL * SDK ( L ) » I NT , 068 

3P( I, 1), P( 1, 2). P< 2, 1 ) . PI 2*2 ) ) 069 

GO TO 168 070 

C SEGMENTS K AND L SHARE UNF POINT (THEY INTERSECT) 071 

80 KG=0 072 

JM=KB 073 

JC =K A 074 

K F - 1 075 

I N f ) = C KR-t. A )* ( KR-LH ) 076 

I F ( IND.NF.OIGO TO H2 077 

JC=KB 07 ft 

KFc-1 079 

JM-K A 060 

KG* 3 08 1 

82 L G* 3 082 

J P s L A 003 

LF*-1 08 A 

IFCLR.EO. JOGO TO H3 085 

J P = L B 086 

L F * l 087 

L G-0 088 

83 SGN=KF*LF 089 

CPSI*( (X( JP| -X( JC) )* ( X ( JM )-X( JC )>♦( Y( JP )-Y( JC ) I* (Y(JH)-YIJC) ) 090 

2*( J ( JP)-Z ( JC I KM l ( JM)-/ ( JC ) ) )/(OK*Ot ) 091 

CALL Z GMM ( .O t OK, .OtOL • AK.COKCK ) « SDK ( K ) • SDK ( L ) » CPS 1 0^2 

2«C( 1« 11*0(2.2) I 093 

on 98 KK = 1 .2 094 

KP=! ABSCKK-KG) 095 

00 9H L L = 1 t ? 096 

LPe I ABS( LL-LG) 097 

PC KP. L P ) =SGN*U( KK. LL ) 098 

98 CONTINUE 099 

GO TO 168 100 

C K=L (SELF REACTION OF SEGMENT K) 101 

120 3 - . 5 102 

1 F ( K A • N 6 • I A ) S=-. 5 103 

CALL 7 GMM ( .0, OK * OK <• ( . 5-S ) * 0K$ ( . 5*S I . AK f COK ( K ) , SOK ( K ) , SDK (K) , 1 . 104 

2*P(l*l)*P(lt?)*P(2*l ) » P ( 2. 2 ) ) 105 

168 CI.J = F1*FJ>P( IS. JS) 106 

I F ( J. GT • NCM ) GO TO 190 107 

C( 1, J)=C( 1 .J)*C1J 108 

I f < i.nf.jici j. n=cc j. n+ru - 109 

GO TO 200 1 10 

190 JG* J- JPP 1 1 1 

C ( I » J G ) = C ( 1 » J G ) -C I J 112 

200 CONTINUE 1 13 

ZH*|.0«.0I H* 

IFCCMM.LE.O. IGO TO 262 115 

CALL ZSURFf AK,GMM f FMCfZS) 116 

ZH=ZS/( 4,*TP*AK ) 117 

DO 260 K=1 , NS 118 

NDK = ND ( K ) 119 

ZSAH*2.*ZH*(0(K ) - S OK ( K )*COK( K M/SDKC K 1**2 120 

00 210 11 = 1* NUK 121 

1 = MU ( K . I l ) 122 

210 IF( !.LF.NCM)C( I. U*L( It II^ZSAM 123 

IFCNPK.EO. i)GO 10 260 124 

Fig. 8. Subroutine IDANT - continued 
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ZOPP=?.*ZH*l SDKCK)-DCK)*COK(K) )/SnK(K)**2 125 

K A = 1 A ( K ) 126 

K8=IB(K) 127 

DO 260 I 1 = 1,N()K 1 2 B 

1 =MD( K, II) 129 

1 F ( I.GT.NCM1GO 10 26 0 130 

F 1*1. . 131 

IF1KB.F0. 12(1) )G0 TO 236 1 32 

I F ( KB. FQ* Ii(I))FI=-l. 133 

15 = 1 1 34 

GO TO 240 135 

236 1FCKA.E0.13( 1 ) >F1=-1. 136 

I S = 2 137 

240 DO 260 J J= 1 ? Nl)K 138 

J=MO(K,JJ) 139 

1 F ( I.QF.JIGO TO 260 140 

F J= l , 141 

I F ( KB • F 0. 1 2 ( J ) ) Gtl TO ?46 142 

1 F ( KB* FO* IK J) ) P J =- 1 • 143 

JS= 1 144 

GO TO 250 145 

246 !F(KA.E0.13( J))FJa-l. 146 

J S= 2 147 

250 1 F < 1S.6Q.JS)CIJ=FI*FJ*ZSAM 148 

1 F C IS.NF«JS)CIj=F I * F J $ Z OPP 145 

IF< J.GT.NCMIGO TO 259 150 

C( I, J)=C( It JKC1 J 161 

C( J. !)=C( J. 1 l+CIJ 162 

GO TO 260 153 

259 JG= J— JPP 164 

CII. JGKCC It JG)-C»J 165 

260 CONT 1 NUF 166 

262 1F(M n.Lri.OlGO TO 300 167 

DO 280 I =1 .NCH 168 

J J A * J A ( 1) 169 

J1=JJA 160 

I 12=121 1 ) 161 

111 *IH ! 1 162 

IF ( I 12. FQ. I B( J 1 ) ) Jl =J1 + NKS 163 

1M I .LE.NPGP1G0 TO 270 164 

J JB= JB ( 1 ) 1*6 

J 2 = J J D 166 

I F f I I2.F0. 1IU J?i 1J?«J?*NRS 167 

C( l« I ) - C ( It 1 )+2LD<Jl)+2LD(J2) 168 

J J J = J J A 169 

DO 268 K = l , 2 1 ?0 

Nf)J = ND( JJJ) 1 ? 1 

DO 266 JJ = 1 t NO J l 7 ? 

J=MD( J J Jt J J 1 I 73 

IF ( J.FO. 1 ) GO TO 266 l 74 

1F( I2( JJ.NE. I121G0 TO 266 1 ? 6 

F 1=1 . I 76 

I F ( K. 60. 2 1G0 TO 264 I 77 

I F I I 1 ( J ) .Me . II l)F 1 = -l. I 78 

CI!fJ>»C<!tJI*M*ZlOCJl) l ?y 

GO TO 266 IPO 

264 I F ( I 3( J ).Nfc. 1 3( 1 ) )H =-l . 1**1 

C( I . J ) =C l I tJ)*F!*ZLD(J2) l e? 

266 CONT I NUt l* 13 

268 J J J= J JB l 04 * 

GO TO 280 166 

270 I F ( IBl JD.IF.NPGP) J1=J1*NRS 166 

Fig. 8. Subroutine IDANT - continued 
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C( It I )=C< If I )+2,*ZlD(Jl ) 16" 

MOJ=ND(JJA) 

00 27 8 J J - 1 t NI)J 

J = W)( JJA, J J ) 1* 

1 F I J.EO. I ) GO TO 278 IV. 

1 F ( 1?( JI.NE.l 12)0(1 TO 278 Is. 

r- 1 = i . i *> 

1 F C IKJI.NF.I II )H=-1. IV- 

CC I • J )*C< I • J > + 2. I*ZLD( J1 ) IV 1 

278 CONTINUE lw ' 

280 CUN 1 INUE 1^ 

300 RETURN 1** 

END 1V‘ 


07 : 

Fig. 8. Subroutine IDANT - continued 


30 


StlAROUT INF IFF l IX I A, !H, INS, I » • I?, 1 3 , MD, N , NO* N S , C OK , C J , D , 00] 

2E PP • F 1 l,FPM # r lH t r, # *CPP,GY 1 , ph, SDK, Th, x, y, 7 ) 00 ? 

COM PL t X l PH, r 1H, CJl , M 1 ,FT?, I PI , FP? C03 

COMPLEX CJ( ) ) .FPP( 1 ) , t 1 T< 1 ) OOA 

0 I MANSION 1 All), Mil 1 ) ,ND( 1 ) ,COK( 1 ) , SDK l)l,D<l),X(ll t Y(l),m) 006 

OlHFNSIfJN l)( 1 ) « I ? I i ) « 1 3 ( 1 ) t MD ( 1 NS , 4 ) 006 

OATA CJI /I . 0* - . S30R8HF-?)/ 00/ 

iHKc.O) /aVXXh 008 

CTHrCnSUHR) 009 

STHsSIMI IHR ) 010 

PHR*. 01 7A633*PH Oil 

CPH=cnsi phr ) oi? 

$PH=SIN(PMR) 013 

00 130 l.-l.N 01 A 

FTT( 1 ) = I • 0, • 0 ) 016 

130 eppms(.o,.o) oi6 

DO 1A0 K*1,NS 017 

K A « I A i K 1 018 

KB= I B ( K ) 0 1 V 

CALL 7 F M X ( K A MIKA), / (KAl,X(KB),Y(KB)*Z(Kh)»OlK) 0?0 

? « CDK ( K ) , SDK ( K ) « C TH, S I H , C PM, SPH, 6 T 1 • E T 2 , FP 1 , EP? ) 0? 1 

NDKeNOIK) 0?? 

DO 1A0 11*1, NO* 0? 3 

1 =MD( K ,11) 0?A 

M - 1 . O? 6 

1F(KA.F0.I?( I 1)00 10 136 0?6 

IF( KH.FO. Ill I ) )M *-l • 0? 7 

£PP( I ) s F P P ( 1 ) ♦ F 1 * F P 1 0? ft 

f T T ( 1 ) = F 1 T( I ) ♦ F I * FT 1 0 ?9 

GO TO 1 A 0 030 

136 IFCKA.f'.'. !3( ! ) IF! — !. 031 

EPPI 1 )=FPPI 1 )+F !*FP? 032 

F T T ( I ) = PI T ( I 1 ♦F I * 6T ? 033 

1 AO CONT 1NUF 03 A 

F PM= ( • 0* • 0 ) 036 

F. T H* ( • 0, • 0 ) 036 

200 00 ?60 I r 1 , N 037 

LTH»FlH^CJf I ) * t T T ( ! ) 038 

?60 FPHrFPH*(M( 1 )*tPP( 1 ) 039 

1 PPcCAMSI F PM ) OAO 

A1T*CABS( FTM) OAl 

GPPsAPPo APP/I 30. *G ) OA? 

GTT*ATT*ATT/C 30,*G) 0A3 

RETURN OA A 

END 0 A 6 

Fig. 9. Subroutine IFFLD. 
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SUBROUT I WE 1 SOR I ( l A « 1 H , ICC* 1 C J • 1 NS « 1 WRITE* 11,12, 1 3 • J A 

• JB, 

001 


?MAX, MIN, v t>, N, NCM, Nu, NP, NS. NkP, NR S , NPGP , NSGP , OC * > C • s C . 

ZC) 

00? 


dimension j mm ? o i » > r. 1 1 1 , / c 1 1 1 • Yt tn,/c(i) 


003 


DIMENSION III 1 >• 1 ?l 1 I* 13(1 ) , JAt ) ) , JBI 1 J 


004 


D 1 MF SION I At 1 ) , IH t 1 ) * r ,0 ( 1 ) ,F*0( 1NS.4 ) 


006 

1 

FORM AT ( NX, • J » • IX, 1 U 1 J ) ' , D i 1 1 8 I J ) • , 1 2X * • K % 


006 


? 1 X , • l A ( K ) 1 , 1 X , • I H f K ) 1 , 7 * , • DC I J ) 1 1 


007 

? 

FORMAT! U, i\ S, 1 I 1 6, 216, 2F1S.6) 


OOfi 

3 

FORMAT I 1 X, 1 I 6, 3F10.S, 1 1 6, 3F 10. b) 


009 

4 

format t i * , m 6, i i i 6, 6 i 6 i 


010 

6 

FORMAT ( Imo ) 


01 1 

6 

TfiRMA T ( SX, • I », 4*, «xc ( 1 ) M^X, • YC < n • ,6X, UCI1) ' • 


012 


2 6X, ' J ' . 4X , 'XC( J I * • NX , »YC 4 J) • * t>X, '/CIJI'I 


013 

y 

FORMAT ( sx, • 1 », 3* . 'JA'.'U, • JB' , 3X, • 1 1 • , 3X , • I ? • , 3X • • 1 3 • 

, 

014 


214*, «K • , ^ x , * J A ' , 3X , • J o 1 , 3 X , • I l • , 3X , * I 2 1 • 3X , • 13* ) 


016 


1 GPPcNPGP* 1 


016 


NP I *NR P-NPGP 


017 

N't- xi SM UP IMF image segments 


0 1 B 


DO IB J- 1 , NR S 


019 


K * J ♦ NR S 


020 


I A l K 1 1 1 A ! J ) 


021 


IF ( I A! J J . GT . NPGP )|A(KI*IA(JMNP! 


022 


I B 1 K ) * 1 H ( J ) 


02 3 

1 B 

I F ( 1 B( J I. GT. NPGP ) !H( K I* IH( J 1«NP1 


02* 

Next SET UP IMF image POINTS 


02 b 


DO 20 I = 1C *P,NRP 


0?6 


J* I 4NP l 


n?7 


xct j » *xc 1 1 ) 


02* 


YCt J)*YCI I 1 


029 

20 

2 C ( J 1 * ~ Z C ( I 1 


030 

NeXT CALCULATE T nr. SEGMENT (.tNGTnS OC ( J ) 


U3 1 


1F( IWRI TF.lt .OIGO TO 22 


03? 


WR ITF t 6, b I 


033 


WRITE (6,1) 


034 

22 

0(1 26 J = 1 * NR S 


036 


K * 1 A { J ) 


036 


L * I B ( J ) 


03 7 


DX*XC(K)-XC(l) 


G3B 


DY*YC( K )-YC( L ) 


039 


D2»?C(K)-2C(U 


040 


OC(v) ) * SOP T ( DX *DX*D Y*DY« DZ * D2 ) 


04 1 


K* J4NRS 


04? 


o r :ki*dc( j) 


04 3 

2b 

»F( JWRITE.GE* 1 )WRlTE(6t2)wt I A ( J ) , l B ( J ) * K , J A ( K > , 1 B ( K ) , 

DC t J ) 

0 A4 


IF ( I WR 1 T 6 • L E . 0 ) <»l) TO 32 


046 


WR I TE ( 6 , b ) 


046 


WR I T F ( 6* 6 I 


04 7 


00 30 1=1* NRP 


04B 


IF ( I .GT.NPGP I GO TO 28 


04 9 


WRITFI6, Ml.xct l ),YCi 1 ),ZCI 1 ) 


060 


GO TO 30 


061 

2 8 

J * I ♦ N P 1 


06? 


WRITE (6, 3) I ,XC( I ) t YC( I ),ZC( lit Jt XC IJ ) ,YC( J) ,ZCi J) 


063 

30 

CONT I NU6 


064 


WRITE! 5) 


066 

CHECK INPUT OATA FOR CONSISTENCE 


066 

3 2 

N* 0 


067 


MIN«100 


068 


MAX. 1 00 


069 


IFINPGP.LE.OIGO TO 40 


060 


no 38 1*1, NPGP 


Obi 


1*0 


062 


Fig. 10. Subroutine ISORT. 
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. { Or’ 'I't.d 

W, 1 ( t !J\ t \ I . X Vivj L> 1*5 t 00 x* 



00 35 J*1,NSGP 

063 


K«( 1A( J)-l )•< IH( jf-l ) 

064 

35 


065 


IF ( L. L T.hax )haa = L 

066 

38 

N=N*?*L-1 

06 7 

40 

iFINRP.LE.NPGPiGO TO 50 

068 


DO 46 I * l GPP « NR P 

06 9 


1*0 

070 


DO 44 Je»,NKS 

071 


K«C I A( J )-!»•( IB( J 1-1 1 

072 

44 

IF <K.F0.0)l«l4 1 

073 


IF(L.IT.M|N»M1N*L 

074 

46 

N«N*2ML-1 ) 

075 

50 

IFfN.lt. 0 .OK. N.GT.ICJIGO TO 500 

076 


IFfMAX.Lt.O • t<R . MlN.Lt. 01 GO TO 500 

077 


IHNPGP.lF.OlGO TO 5 H 

078 

set 

UP TMF MODES Al IMF GROUND PLANE THAT Will NOT HAVE i MA Gt S 

07V 


DO 56 I * 1 1 NPGP 

050 


J*0 

OF 1 

5? 

J*J4 1 

J82 


1 AJa 1 A( J i 

083 


1DJ«1B( J) 

084 


KK*f 1 A J — 1 )*( IHJ-I 1 

085 


1FI J.eO.NSGPIGfl TCJ 5 4 

086 


IF(KK.NE.O)GO 10 5 2 

08 7 

54 

JM 1) *j 

088 


JB( 1 )«J4NRS 

08V 


i?( n-i 

ovo 


III 1 ) * ! ft J 

OV 1 


iff ibj.eo. n 1 1 ( f >«r aj 

092 

56 

13(101111) ♦NPI 

093 

56 

I • nPGP 

094 


NrNPGP 

095 


NCM*NPf,P 

096 


JPPrO 

097 


IMNRS. FO.NPGPiGf) TO 75 

098 

SET 

UP THt -E ST OF 1 H t REAL MODES 

099 


DO 65 K * 1 » NR P 

100 


NJKsO 

101 


DO 60 J-i.NRS 

102 


1ND*( I A( J) -K )o ( I R 1 J)-K ) 

103 


I F < I NO. NE. 01 GO TU 60 

104 


N JK r N JR ♦ 1 

105 


JSPI NJK ) t J 

106 

60 

CONT 1 NUE 

107 


MOO c NJK* 1 

108 


IF (MOO.LF.OIGO TO 46 

109 


DO 6? I MO c 1 • MOO 

110 


l*Ui 

1 l 1 


1 PI). | MO. 1 

112 


JA ! * JSP< I Ml) ) 

113 


J A ( 1 ) s J A ! 

1 14 


JB I r JSP( I PD ) 

115 


JB( 1 )* JBI 

116 


11(1 )«|A( JA1 ) 

117 


V F C IA( JAI ) • E 0. K ) 1 1 ( l ) * I B ( J A I » 

116 


i?( n-K 

119 


I3( I )*l A( JBI 1 

120 

6 2 

IF ( 1 A< JBI ) .EO.K ) I 3( I )*10( JB!) 

121 

65 

CONT I NUE 

122 


NCH-- I 

123 


JPP.NCM-NPGP 

124 


Fig. 10. Subroutine ISORT - continued 
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SET UP IMF IMAGE MODES 


125 

DO 70 1*1 c»p° , ncm 


126 

K * 1 ♦ J P P 


127 

JA( K ) * J A I I MNRS 


12R 

JBSIUeJHIIMNMS 


129 

I I A = 1 1 ( I ) 


1 30 

i ib = i?( n 


131 

1 1 C a 1 3 ( n 


132 

1 1 ( K ) * I 1 A 


133 

IF ( I l A.GT . NPGP ) IK K 1 * I I A + NP.I 


13V 

I ? ( K ) » 1 IK 


135 

IK IIR.GT .NPGP ) 1 ? ( K 1 » I 1 B*NP 1 


1 36 

I 3 ( K 1 - I I C 


137 

70 IF ( 1 IC.G1 .NPGP ) I 3< K 1 = I lONP I 


138 

Ns?»NCM-NPGP 


139 

75 M A X * 0 


1V0 

M I N*. 100 


1M 

NO ( J ) « NUMBER 0* DIPOLE MODES 

SHARING SEGMENT J 

IV? 

HOI J » K ) * LIST of UK'ULES SHAKING SEGMl NT J 

1 V 3 

DO 100 J 3 1 * N S 


1 <*<» 

DO 80 * * 1 , V 


1 V5 

80 MD ( J* X ) =0 


1 V6 

K*0 


1 V 7 

DO VO I *1 ,N 


1 V6 

J A 1 s J A I 1 } 


1 V 9 

JOI =JH( I I 


150 

l «I JA!-J)M JflI-JI 


151 

IFU.NF.OIGO 10 VO 


152 

K*K* 1 


153 

MD( J, K ) * 1 


15V 

vu Coni inue 


155 

NO ( J I 


156 

If 1 K.Gl . MAX )MAX *K 


157 

100 IF ( K.LT.MINIMINsK 


158 

IF ( lMR1If.LF.0IGn 10 500 


159 

MR 1 T f ( N* 7 1 


160 

DO MO l-I.NCM 


161 

IF ( I.GT.NPGPIGU ID 108 


16? 

WRITE (6tM I * J A I I ), JH( I )• I 1( 1 

1,12(11,13(1) 

163 

GO TO 110 


16V 

108 K«tOPP 


165 

WR1TF(6,A)|, JAC I ) • J B ( 1 ) * ! 1 ( l 

)• 121 I), I3( 1 ),K, JA(K) v JH(K) , 

166 

?llCKltl?IR)*I3(K) 


167 

1 1 0 COST INUE 


168 

WRITE <6,51 


169 

500 RETURN 


170 

END 


171 


Fig. 10. Subroutine ISORT - continued 
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SUBROUTINE PDI SSI I a. 18. INS. M • 1 ?• 1 3,MD.ND,K'l D.NS • C J . C*m, 0. CDK , 001 

?SD< • D ! SS • 2H f 2 1 0 ) 00? 

COMPLEX CJ< 1).Zh,CJA,CJB./LA,ZiB.7LD< 11 003 

dimension :o < m ksom n.ou > . i i i i » . i?u i , i 3111 . i ai n . ibi n .Nom ooa 

DIMENSION m:j(|nS,a) oob 

RM.RFAL 1 ZM) DON 

D I S S* . 0 007 

on ioo ks i ♦ n s oort 

IMChm.ll.o. K.ti TO 10 . OCW 

F A * ? . • KM» ( ll| K 1-SlKI K|*CI)K|K| l/SOKIK )*•? 010 

FB«A.*RM*| SDK! KH)K)0CIXUJ l/SOKI < )••? 01 1 

10 KA>|A|<) 01? 

Kf)tlH(Kl 013 

CJAt|.0,.0l 01'. 

CJO«(.0,.0) 01b 

NDKtNDIK) Olfc 

00 AO II » l.NOK 017 

1 «*mi) ( k . 1 1 ) oi a 

f i ■ l . 01 V 

IF I KB. ( 0. 1 ?l 11)00 TO 36 0?0 

Iff KH.FQ. 1 1 111 )F I «-l. 0?1 

CJA«C JAif I *CJ( 1 ) 0?? 

oo in ao o ? 4 

36 IF I KA.FO. 1 3C 1 ) >H *-l. 0?A 

C JB^C.IB^F 1 *C Jl I > 0 ?b 

AO CONTINUE 0? 6 

IF I MID. IF . 0)00 10 SO 0? 7 

AJA-CAHSIC JA )•♦? 0? B 

BJB-CAHSICJtt)**? 0?V 

KK»K«NS 030 

RtA*RFAll 7101* 11 C31 

RLB^RF AL ( ZL1M KK ) 1 03? 

01 SS*DI SS^AJA*«L A4HJfl»RLB 033 

SO IFCCMM.LF.O. IOO in IOO 03*. 

r)!SS«DISS*FAc{CAftSICJAM*?4CAHS<CJH)**?l 0 3b 

?#Ffl* ( at All c JA )*Ht luc JR) ♦A I MAC.I C JA )*A 1 MAGIC J8) I 036 

100 CTINT IMit 037 

RETURN 03B 

END 03V 

Fig. 11. Subroutine POISS. 
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SUBROUTINE ZEE ( XA, YA, 2 A, XB, YB, 28, D 001 

2, CKO, SKD, CTh, S TH, CRH, SPH, FT 1 , FT 2, EP1 • EP2 ) 002 

COMPLEX EJA f EJB,tPl,EP2,ESl,ES2,ETl,E12 003 

CA = ( XB-XA )/l) 004 

CB = < YH-YA ) /I) 00b 

CG= ( 2 B-Z A ) /D 006 

G=(CA<‘CPH-»-CB‘ > SPH)*STH4CG*CTH 00 7 

GK= 1 • -G*G 008 

E T 1 = ( • 0 , . 0 ) 009 

FT2=(.0,.0) 010 

EP 1 = ( .0, .0) Oil 

E P2= ( • 0, . 0 ) 012 

] F ( OK . L T , . 00 l ) GO TO 200 013 

A = XA*STH*CPH + YA*STh<'SPHWA*CTH 014 

B = XB#S1 H*CPH + YH#S I h*SPH4 7 b*CTM 01b 

E JA=CMPLX (COS( A) , S1N( A) ) 016 

EJB=CMPIX(C0S(B),51N(B)) 017 

SGD=SIN<G*D) 018 

CGO=COS ( G**0 ) 019 

E S 1 =30. JA*CMPL X ( SGO-G*SKOf CKD-GGD ) /GK/5KD 020 

ES2 = 30.*E JH^CMPIX < G-SKU-SGO,CKD-CGO) /GK/SKD 02 1 

T*( CA*CPH+CB«SPH ) *CT M-CG*S TH 022 

P=-CA»SPH4CH«CPH 023 

ET1=T*ES1 024 

ET2 = T* E S2 02b 

EP1=P*ES1 026 

EP2=P*hS2 027 

200 RETURN 028 

END 029 


Fig. 12. Subroutine ZFF 


SUBROUT I N6 7GMM( SI , S2, T 1 , T?, U, CGDS, SGOl , SG02.CPS1 , PI 1 , P12, P21 , P22) 001 

DOtllili PRECISION Rl«K2.0P0 t S!S«TSlf lS2tSTl*ST2«r.U v BOtCPSS»SPSltSK 002 

2t TL 1 * Tl.2» TOl ♦ TD2 • SOI , OPS I , DO, ZO 003 

COMPLEX E<2,2),F<2,2),6AM,P11, Pl?tP21 f P22 004 

COMPLEX EB,EC*h*K,El ,FK1 ,FGZ1,FS1, fc S? , E T 1 , H 1 2 , F XP A , E X PB 003 

COM PL C x EG 2 < 2 . 2 ) . G.-U 2 ) • GP ( 2 ) 006 

DATA ETA, GAM, P I /376. 727, ( . 0, 1. ) , 3. 14139/ 007 

DSO=D*D OOB 

SGOS- SGOl % 009 

IF< S2.LT .SI 1 SGDS-- SGD 1 010 

SGOT = SGD2 Oil 

IF< T2.LT.T l 1 SGOT *- SG02 012 

1 F ( ABStfPS! ).GT.0.997)G0 TO 110 013 

ESI=CFXP(GAM*S1 1 014 

ES?=CfXP( GAM$S2 ) 013 

ET 1 =CFXP(GAM*1 I ) 016 

ET2=CEXP( GAM« T 2 ) 017 

00 = 0 018 

OPSI=CPSI 019 

TO 1 =1 1 020 

T 02 = T 2 021 

CPSS=DPSI*L>PS1 022 

CO=DQ/OSQRT ( 1.00-CPSS) 023 

C = C D 024 

BD=C0*0PS1 023 

B = BD 026 

EB=CEXP(GAMcCMPL X( .0, B) ) 027 

EC=CEXP( GAM^CMPLX ( .0,C M 028 

DO 10 K = 1 ♦ 2 029 

00 10 L = 1 , 2 030 

10 t(K|L):(>Oi*0) 031 

TS1=TD1*T01 032 

TS2=TD2*TD2 033 

DPO=DO*DD 034 

SI=S1 033 

00 100 1=1,2 036 

F I = ( — 1 ) ♦* I 037 

S01=S1 038 

SIS=SOI*SOI 039 

ST 1=2.* SOI * TOl- OPS I 040 

ST? = ?.< J SDI*TQ2 <t DPSI 041 

R1=0S0RT (DPO+SlS+TSl-STl , 042 

R2=DS0RT< OPO+S IS+TS2-ST2 ) 043 

EK*EB 044 

00 30 K= 1 , 2 043 

FK«( -1 ) $ * K 046 

SK=FK*SOI 047 

EL=EC 048 

00 40 L = 1 , 2 049 

EL = ( -1 ) 0 * L 030 

EKL=FK*FL 031 

XX=FK*BO+FL*CD 032 

TL1=FL*T01 033 

TL2=EL^T02 034 

RR1 =R 1 *SK4TL 1 033 

RR2=R2*SK*TL2 036 

CALL EXPJ( GAK*CMPLX< RR 1 , -XX ) , GAM*CMPL X I RK2 , -XX ) , EXP A ) 037 

CALL 6XPJ(GAMoCMPLX(RRl, XX) ,GAM«CMPLX(RK2,XX) ,EXPB> 038 

E(K,U=E(K,U*Fl* ( 6XPA*EKL^FXPB/EKL ) 039 

40 FL=1./EC 060 

30 E K= 1 . /E 8 061 

ZO=SOI*OPSI 062 


Fig. 13. Subroutine ZGMH. 
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zc = zo 

EG Z l=CFXP(GAM*ZO> 

RR 1 = R 1 ♦ ZU- T D 1 
RR? = k? + /A)-TOZ 

CALL FXPJ(GAMokR l , GAMCRR2, EXf>B) 

RR 1 =R 1 -?.[)♦ mi 
RR?=R2-Zn*T02 

CALL EXPj( GAM*RK1 f GAK*KR2, tXPA ) 

F ( 1 • 1 ) = 2. * SGOS* ( .Ot 1* )*EXPA/EGZ 1 
F ( t ,?)=?.« SGDS* (.0,1. )*EXPI)#6GZ I 
100 Sl=52 

CST--ET A/( 16. *P i*SGOS*SGDT ) 

P 1 1 =CST * ( ( F( l, 1 nE(2,2)*ES?-E( l»^)/ES2)^FT? 

A 4( -F(l • 2)-t<2, 1)*ES24E( 1 . n/ES?)/fc12> 

P12=CST*( (-FU, i )-E ( 2,2 )*FS?*E( 1, 21/ES2 )*ET 1 
B 4( F( 1 , 2 ) 4H 2 , 1 ) *ES2-E ( 1. 1 )/ES2)/Fl 1 ) 

P21sCST* ( (~JF(2, 1 ) -F( 2, 2 )<ESl4E( I. 2 1/tSl )< fcT? 

C 4( F(2,2)4E(2, i ) * E S 1 - E ( 1, 1 ) / E S 1 1/ET2) 

P22-CS7 * ( ( F(-2, 11 4f( 2, 2 (♦FSl-Ef 1 • 2 1/tSl 1*F. T1 
0 * (-F ( 2, 2)-fc(2, l )*ES14E( 1, 1 1/ES1 l/ETl ) 

RETURN 

110 1MCPS1.LT. 0. )Gl) TO 120 
T A = T 1 
1 B~T 2 
GO 10 130 
120 1 As-1 1 
1 R = - T ? 

SGD1 =-SG()T 
130 S1=S1 

DO ISO 1*1,2 
T J = T A 

DO 1A0 J = 1 , 2 
Z I J=T J-SI 

R=S0RT(^S0>7 J J^ZIJ) 

W^R+Z 1 J 

1F(ZIJ.LT.0. )W = l)SO/(R-Z 1 J 1 
V- R-Z l J 

1 F I ZI J.GT.O. ) V = DSO/ ( R 4 Z I J ) 

IF (J. ED. 1 )V1=V 
IF( J.EO. 1)W1=W 
EGZ ( I , J ) = CE X P ( Gam*Z I J) 

1 AO TJ-TH 

CALL EXPJ ( GAM* VI , GAM* V, GP( 1 ) ) 

CALL EXPJ( GAM‘W1 ,GAM«W f GM( 1 ) ) 

150 S J =S2 

CST=FTA/(R.*P 1*SG0SoSGDT ) 

PI l =CSTM GM( 2 ) 4fcGZ( 2. 2) 4GPI 2 1/EGZ ( ?• 2 1 
2-CGCS* ( GM ( 1)«EG2( 1 , 2 14CPC 11/FGZ ( 1, 2 1 > 1 
P12 = CS1M -GM(2)*EGZ (2,1 1 -GP( 2) /FGZ (2,11 
2 4 0005“ ( GM( 1 ) - FGZ ( 1 , 1 ) 4(,H ( 1) /EGZ ( 1 , 1 ) ) 1 
F21 -CST*( GMI 1 ) ► 0 7 ( 1 , 2) 4GP( 1 ) /MiZ ( 1,2) 

2-CGDS* ( GM( ? { 2, 2 ) 4GP( 2 1/fcGZ ( 2, 2 1 1 1 

P?2=CST“( -GM { I )« EGZ ( 1 , 1 )-GP( 1 1/EGZ ( i, 1 1 
2*CGDS* ( GM ( 2 ) * EGZ ( 2 * 1 ) *GP ( 2 ) / EGZ l 2 , 1 1 1 ) 

RETURN 

END 

Fig. 13. Subroutine ZGMM - continued 
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06 A 

066 

066 

G67 

068 

069 

070 

071 

07? 

073 

07A 

075 

076 

07 7 

078 

07 9 

oeo 

08 ^ 

082 

08 3 

08 A 

085 

086 

08 7 

088 

089 

090 

09 1 

09? 

093 

09 A 

095 

096 

09 7 

098 

099 
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103 

10A 

105 
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107 

10b 

109 

no 
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SUBROUT IMF ZGS( XAA A, Z A, XB, YB, ZB, XI ,Y1 , Z 1 *X2, Y2* Z2, AK, 001 

2DSrCDS,SI)S,()T*SUf* I NT , P 1 1 , P12, P2 1 , P22 ) 002 

COMPLEX CST* FJ1. E J2, EJA, EJH, FR 1 , EK2, ET 1 , ET 2, PI 1 , P12, P21 *P22, GAM 003 

COMPLEX SGUS,SGI)T OOA 

DATA fclA,GAM f Pl/376. 727, <.0, 1. ), 3. 1A159/ 006 

C A = ( X 2 -X 1 ) /OT 006 

CB=( Y2-Y1 )/Dl 007 

CG=(Z?-71)/0T 008 

CAS- ( XB-X A ) /OS 009 

CBS=( Y B - Y A )/D$ 010 

CGS=( ZB-ZAI/DS 01 1 

CC=CA*0AS*CB*CBS4CG*C(»S 01? 

IF ( ABS( CC ) . GT. 0. 99 . ) ,U TO 200 013 

20 SZ=MX1-XA)*CAS+(Y1-YA)*CBSMZ1-ZA)*CGS 0 1 A 

1 F ( IN7.E0.0JGU TO iOO 0 1 6 

COOS -COS 016 

SCO S- CM PL X ( . 0, SOS ) 017 

SGDT=CMPLX( .0,SDT ) 01B 

I N S - 2 $ ( 1NT/2) 019 

IF ( INS.LT.2) INS-2 020 

I P= I NS+ 1 021 

OELT=OT/INS 022 

T = . O O? 3 

OSZ = CC*()EL T 02 A 

Pll=I.O t .O) 026 

Pi?=(.0,.0) 026 

P? 1 = ( • O, . 0 ) 027 

P?2 = ( . Of . 0 ) 028 

AKS=AK*AK 029 

SGN = - 1 . 030 

DO 100 IN=1 , IP 031 

ZZ1=SZ 032 

Z Z 2 = SZ-OS 033 

XXZ-X1+T*CA-XA-SZ«CAS 03A 

YYZ cY 1 ♦T*CR-YA-SZ*CBS 036 

211 = Z 1 ♦T*CG-ZA-SZ*CGS 036 

RS = XXZ$o2*YYZ<- rt 2 + ZZ Z**2 037 

R? =SORT < KS + ZZ 1**2 1 038 

EJA=CMPLX(C0S(K1 )t-SIN(Pl ) ) 039 

EJ1«EJA/R1 OAO 

R2 = SOR T (RS*?Z2**2) 0A1 

EJB-CMPLXIC0S(R2),-S1N(R2) ) OA? 

E J2~E Jft/R 2 OA 3 

ER1 =E JAo SOUS ♦ Z Z l#fJl*CGDS-ZZ2*EJ2 OAA 

ER?«-E JB*SG0S*ZZ2*EJ2»CGDS-ZZ 1*EJ1 0A5 

F AC - • 0 0A6 

I F ( KS.GT. AKS )F AC = <CA*XXZ+CB* YYZ *»CG*ZZZ )/RS 0A7 

ET 1 =CC*(EJ2-E J1*CGUS ) *EAC*FRl 0A8 

ET2*CC* : ( EJ1-EJ2*CGUS)+FAC*ER2 0A9 

C = 3 • ♦ SGN 060 

I F { 1 N * E 0 . 1 .OR. IN.EO. I P ) C = 1 . 061 

Cl =C*S 1NI0T-T ) 062 

C2 = C° S 1 N ( T) 063 

P)l=Pll*FTl*Cl 06 A 

P12=P]?+ET1*C2 066 

P? I = P2 l ♦ F T ?*C 1 066 

P22-P22 + E T2*C2 067 

T=T4UELT 068 

SZ=SZ*DSZ 069 

100 SON = -SON 060 

CS1=-(.0, l.)*ETA»OtLT/( 1 2. * P I * SOUS* SGOT ) 061 

P 1 1 =CST <* P 1 1 062 

Fig. 14. Subroutine ZGS. 
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P12=CST*P12 063 

P21=CST*P21 064 

P22=CS1 *P?2 066 

RFTUKN 066 

200 S? 1=( X1-XA)*CAS>( V? -YA)*CHS*( 21-ZA4*CGS 067 

RH1 =S0RT ( < X 1-XA-SZ l*CAS )**2*< Y 1-YA-S2 1*CBS>**2*( Z 1-2A-SZ 1*CGS ) **2> 068 

S72 = SZl+l)T*CC 069 

RH2*= SOR T { ( X2-XA-S?2*CAS)**2+< Y2-YA-SZ2*CBS)**2-M Z 2-2 A-SZ 2*CGS ) **2 ) 070 

00K=( RH1+RH2 )/2. 071 

JFfDOK.GT. 20. OAK .AND. INT.GT.OJGO TO 20 072 

IF ( ODK. LT. AK )Dt)K = AK 07 3 

CALL 2GMM( .0,0S»SZ l , SZ 2, DDK , CDS, SOS, SOT , 1 . , P )1 , P 1 2 , P2 1 » P22 ) 074 

RFTUKN 076 

300 SS=S0RT( l.-CC*CC) 076 

CAD=(CGS*CB-CHS*CG) /SS 077 

C B D = ( C A S * C G-C G S * C A ) / S S 07 8 

CGD*(CflS*CA-CAS*CB)/SS 079 

DK = I Xl-XA )*CA1>4 l Yl-YA )*CB0*( 2 1-ZA)*CGD 080 

DK=AHS(OK) 081 

1FI0K.LT. AK)0*=AK 082 

XZ=XA*SZ*CAS 083 

YZ=YA4-S2*CBS 08 4 

Z7«*ZA*S2*CGS 086 

XPJ =X 1 -DK*CAD 086 

YP1 =Y 1-0K*CBD 087 

2P1 =2 1 -DK*CGO 088 

CAP=C8S*CG0-CGS*CB0 089 

C H P =CG S * C A D-C A S <*CG D 090 

CGP=CAS*Cm)-CHS*CAU 091 

PI *C AP*( XP 1-XZ ) 4CBPM YP1-Y2 ) *CGP*< 2 P 1-22 1 092 

T1-P1/SS 

SI =1 l^CC-SZ 094 

CALI 2GHM ( S 1 * S 1 + DS» T l.Tl*DT,DK t Cl)S*$DS* SOT , CC , P 1 l , P 1 2 , P2'l • P22 ) 096 

RFTUKN 096 

f NO 097 


Fig. 14. Subroutine ZGS - continued 
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SHRROU T 1 N8 ZSURM AK , CMM , F^C • l S ) 

COMPt F X BFS.HrS 1,/S 

DMA PT A , SOI , TP / s?/>. IV 121 , 1 . 4 14? 1 356, t> . 2H31 hb3/ 

SOSwi * ) . f 6* SDK T ( (. M M/ 1 P/FKC/H. 8S43 3 ) 

X = AK v SDSv.( 

If ( X.r.T.H, )C>n 7 : J SO 
7 = X / R • 

T2&TM 

T4rT?*T? 

BE k = ( ( I ( ( {-.90)1 -5 • T 4 ♦ . 12?bS?P~2) *7 4-. 03349609 )*T4* 2.64 19 14 )*T4 

2-3 ? . *63*. **6 J * T 4* l | 3 . 7 7 778 ) ' T 4-&4 . ) ^ U- ♦ 1 . 

fit ] ... ( ( ( ( ( ( . i i 34'. P- -i • T 4-. 0) l 0 3 667 ) * 7 4 ♦ . b? 18 561 S i * I 4- i 0. /6 r >b ) * T 4 

2 « 7 2 . R 1 7 7 7 7 ) rt T 4 - 1 1 * . 7 7 7 7 H ) <• I 4 <• 1 ' . ) e 5 2 
h . ; k$T 2«(U(H . • >41 - •;!•.«. 4b9b 7F-3M i 4-.02609253 ) *• • ••* • ••*' t 4 784V ) 

2*74-6. 06 H 1 4 H 1 )v 144 14.222222 ) '*14-4 . ) 

UP I P - X* ( { ( ( ( ( . 4 009F-4M4-. - 7 9 3 K 6 f - 2 ) * 7 4* • 14677204 ) M 4 - ? . 31 167b) ) * 
PDh 11.3/7 ' M ) I 4- 1 o. 56666 7 J * 7 4 * . b ) 

Hf S-CMPLX ' 11 fcRi HP 1 ) 

R { S 1 ~ . / 0 7 1 O * * C M v t X ( b b R P - B l : I P t 8 b R P ♦ R C 1 P ) 

GO TO 100 

SO XP=.7071 06 8 I « X 
Xl=l./X 

K r { ( - . 04 4 9 M*S^v 1 ♦ . * 9 0 6 ?5 F — 2 1 * X 1 ♦ . 0 H P $ B R 3 b ) * X 1 ♦ 1 . 

1 -( ( 04M)3bb9< <1-. 06?b )*X 1-. OHH3HH3b ) K X 1-. 39?6v90 7 *XP 

fU S r f cC MPL X ( Cns I 1 ) , s I N ( 7 ) ) 

F = ( ( . 1 1 2**023 1 C / 1 * . 0? b 1 b6?b J * X 1 - * ?M> ] 6b05 ) * > 1 * 1 . 

7 - ( ( . 1 1 600 / / :: * 1 1 *- 7b ) *-X 1 ♦ . 2 6 S 1 650b ) *X 1 ♦ 1 . 1 7809 /?♦ XP 

BPS 1 *F«CMPL X ( COS ( I ) , S INI 1 ) ) 

100 Z S -- -C.MPL X I l . * - l . )•»■ T A*B£S/BbS 1 / SOT /SOSwb 
Rf TURN 
END 


001 

002 

003 

004 

005 

006 
007 
006 
009 
CIO 
01 1 
012 

013 

014 
01b 
016 
017 
Cl « 

019 

020 
02 l 
022 

023 

024 

025 

026 
02 7 
028 
029 
0 30 
031 


Fig. 15. Subroutine ZSURF. 
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ABSTRACT 


A computer program is presented for a thin-wire antenna over 
a perfect ground plane. The analysis is performed in the frequency 
domain, and the exterior medium is free space. The antenna may have 
finite conductivity and lumped loads. The output data includes the 
current distribution, impedance, radiation efficiency and gain. The 
program uses sinusoidal bases and Galerkin's method. 
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IV 


I. INTRODUCTION 


Reference 1 presents a computer program and reference 2 
presents the theory for thin-wire antennas and scatterers in a homo- 
geneous conducting medium. The present program differs from 
reference 1 only in the following details: 

a. The exterior medium is free space. 

b. The antenna is situated over a perfect ground plane. 

c. The wires have no dielectric sleeves. 

d. The frequency is real. 

e. Scattering problems are not considered. 

To avoid unnecessary duplication, it is assumed the reader is 
familiar with the program in reference 1. 

The program handles antennas constructed of straight wire 
segments. One or more segments may connect to the ground plane, or 
the antenna may be situated some distance away from the ground 
plane. No segment has both endpoints on the ground plane. This 
program can readily be modified to handle more general situations 
as in reference 1. The program uses the delta-gap model for the 
generators. 

The method of images is employed to reduce the problem to a 
thin-wire structure in free space. The user sets up the geometry 
of the real wire configuration, and the program automatically sets 
up the image. I-f unlimited storage were available, one might set up 
a large mutual-impedance matrix for the wire antenna and its image 
in free space. Instead, this program takes advantage of the ground- 
plane symmetry and sets up the compressed matrix C(I,J). Only the 
currents on the real segments are treated as independent unknown 
quantities, and the image currents are dependent. All the currents, 
however, are influenced by the mutual couplings among all the seg- 
ments including the images. In taking advantage of the mirror 
symmet ry, we lose the advantage of having a symmetric matrix. 
However, the net gain in computational speed ana storage is 
substanti al . 

In practice, many wire antennas operate over a ground plane 
with finite conductivity and finite extent. In many cases, however, 
one may substitute a perfectly conducting ground plane of infinite 
extent without unduly disturbing the antenna current distribution or 
impedance. After these quantities have been calculated, one may 
then take into account the finite ground plane in calculating the 
efficiency and patterns. The present program, however, assumes an 
idealized ground plane throughout. 

The remaining sections present the computer program with 
enough explanation to enable an experienced engineer to use it. 




II. THE INPUT DATA AND SUBROUTINE IWIRE 


Fig. 1 is a fortran listing of subroutine IWIRE. This 
subroutine is used to set up the input data. The following data 
must be read or programmed in IWIRE: 


AL 

wire radius a/ A 



CMM 

wire conductivity in 

negamhos/ni 


DPH 

increment in far- fie' 

Id angle <* in 

degrees 

FMC 

frequency in MHz 

TH 

elevation angle e in 

degrees for 

far-field pattern 


To define the shape and size of the wire antenna, the input 
data includes the coordinates XC(I), YC(I), and ZC(I) of the wire 
endpoints, terminals and other current sampling stations along the 
wire axis. The unit of length is selected by the user, and SCALE 
is the conversion factor such that XC(I)*SCALE gives the coordinate 
of point I in meters. NPGP denotes the number of points on the 
ground plane, and NRP is the number of real points including those 
on the ground plane, coordinates are supplied only for the real 
points. The ground plane coincides with the xy plane. 

NRS denotes the number of real segments, and NSGP is the 
number of real segments naving an endpoint on the ground plane. 

For each real segment J, the input data specifies the endpoints 
IA(J) and IB ( J ) . In assigning index numbers to the segments, the 

lowest numbers must refer to those having an endpoint on the ground 
plane. In assigning index numbers to the points, the lowest numbers 
must refer to those on the ground plane. 

Set IWRCJ = 1 to obtain a writeout of the .rtenna currents; 
otnerwise IWRCJ = -1. Set I WRITE = 1 to obtain a lteout of the 
antenna geometry; otherwise IWR5TE - -1. If INT = 0, the rigorous 
closed-form expressions will be used for the mutual impedance of 
sinusoidal dipoles and the calculations will tend to be slow but 
accurate. If INT is a positive integer, Simpson's rule will be 
used for the mutual -impedance calculations. The closed-form ex- 
pressions are always used automatically for the most critical 
impedances. We usually use INT = 0 for multi -turn lo^p antennas 
with closely-spaced turns, and INT = 4 for general purpose. 

Simpson's rule uses INT integration intervals. Thus the accuracy 
and the execution tine tend to increase with larger values of INT. 

Set NLD = 0 if there are no lumped loads; otherwise NLD = 1. 
ZLD(J) denotes the impedance (in ohms) .nserted in segment J at 
endpoint !A(J). A lumped 1 oad at endpoint I B { J ) is denoted by 
ZLD(J + NRS). The user sets up on.y the rea 1 generators and lumped 
loads, and the program takes care of the images. 


2 


StlHMHJT INF 1 w 1 *' * ( \ A, 1 H, ir.v, 1 \S, ; M , 1 .^C J, ! W« II F .Mfi.N?. NS.NKP, 001 

?np s , Npr,p . n * (.»• , /. i , c wm, tvu , f * C . m M 1 # i h, VO, xc » VC • tC • 1 l i ) 00? 

0 1 Of r.'S I ON 1/1(1 }JH ! ),*U l j.tci U./CI 1) 003 

fx VG ( 1 I , / LIM 1 ) 006 

6 H ■ ^nu.'l* ' ♦ l 6 f bx f 'I.MH, • INS- *, J6 t ‘>x, »NS = *t 16) 00b 

b H)kKM i 1*0 006 

C I A | J ) A f:t) Hm .« ) A** INOPCjINTS fit SK.VrM J 007 

C X C ( I ) • Y i f K | v / ( ( 1 I AKf .M t S »»F PCI 1 I i ~1TH AKhllPA^Y IlMTS O'tH 

f. *4*P - \ iH Kl Al POINT',* 1 \C | • ii) J \G Im’.SI « 'i luf GkIHJM) f'LAOfc 00 V 

c OKS a MJMnfv III K | A I st*r» M lNTS 010 

C AllMil a m ht-k UP -.'IMS IN fi*r (KMlNt) »M AM Oil 

c f.sr.p a WiMHtk (It Pt/.J bt-CStMS .. n M fc M)Pl> l NT l.N IK(IUM) PI ANt 01? 

c Nil) =• Nt *hF (if LOADS 013 

Dfl 10 J a 1 • 1 NS 016 

VC( .1 I M .0, . 0) ois 

10 /li e II I .O ) 016 

C SIT UP lHl KPAl ID M «■- A TOWS AND kf-Ai LUMPfO luADS 017 

jgn- 6 oi h 

Vr.(Jf.Msll.iO, I 01 D 

N L 0-0 0?0 

1 *4l t 6 0? 1 

IWHIU- 1 0?? 

IWKCJ'l 0? 3 

AL-.O001 0?6 

CKMa J . 0?b 

l>PH-?0. 0? 6 

fMCaVh. 0?7 

SCAIM1. t?K 

0?9 

NSGPr? 030 

fur.u-2 C31 

NKS*6 03? 

Nkpr6 033 

NS = ?*NPS 036 

NP r P-MMP 030 

Wk I 1 < ( 6 . -• > I V . V , 1 NV.NS 036 

IHNS.GT.INS . »l*. NP.GT.1NP100 U) 600 03.7 

f)fl ?l) 1 = 1. NK »» 036 

XC( 1 1 a.O 03V 

YCllle.0 060 

?0 /C (1 I e . 0 °*»1 

C NFX1 Ml IIP TmF PbAL POINTS 06? 

XC C 1 > = 1 . 04,3 

XC(3)=1. 06. 

?C(3) .b 06b 

ICi6)*.S 066 

XC,(S)*l./07 06/ 

?C(bl-l.?07 068 

c Nf X I Ml Itp Tm| k f A l SeOHfcNIS 06V 

1 A I 1 J s 1 Ob 0 

lhll 1*3 • «bl 

i A C ? 1 - ? Ob? 

I8(?)r4 Ob 3 

I A I 3 1-3 C)t,A 

I #i ( 31*6 Cbb 

IAIM- ; Ob6 

1 I* ( 6 1 r S Ob. 

600 PI UlPN Ob 6 

l NO ObV 


Fig. 1. Subroutine I W I RE . 
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The generator location is defined by JGN . The numbering 
system for the generators is the same as for lumped loads. If the 
generator is to be inserted in segment J at endpoint I A( J) * the 
generator index JGN is the same as the segment index J. To insert 
a generator at endpoint IB(J), set JGN = J + NRS. VG(JGN) denotes 
the complex voltage of the generator. The reference direction for 
these voltages is from IA(J) toward I B ( J i . If the ar.tenna is fed 
with several generators, delete JGN and merely input tr.e generator 
voltages VG. 

f’P and NS denote the number of points .and segments, re- 
spectively, for the complete system (antenna arid image) in free spcce. 


III. THE MAIN COMPUTER PROGRAM 

The main computer program is listed in Fig. 2. Tnis program 
calls subroutine IWIRE for the input data. Then it calls subroutine 
ISORT to generate and store the data for the image points and image 
segments and the length DC(J) of each segment, then ISORT generates 
a list of sinusoidal dipole modes for the complete system (antenna 
and image) in free space. Dipole mode 1 has segments JA(I) and 
JB(I), terminals at point 12(1), and endpoints 11(1) and 12(1). 

This subroutine also generates the following information: 

ND(J) number of dipole nodes sharing segment J 
MD( J,K) list of dipoles sharing segment J 
NCM size cf the compressed matrix 
N number of dipole nodes on the complete system 

The following quantities must be specified in the main program: 

ICC dimension for the compressed matrix C(i,J) 

ICJ dimension related to number of dipole modes N 

INP dimension related to number of points NP 

INS dimension related to number of segments NS 

In Fig. 2, all quantities having the sare dimensions are dinens * d 
in the sane or adjacent statements. The numerical values a *; s i • r • 
to ICC, ICJ, INP and INS must agree with the dimensions actually 
a* served for the corrospondi ng quantities in tne COMPLEX and 
DIMENSION statements. ICC, ICJ, INP and INS must be at leas: cr 
large as NCM, N, NP and NS, respectively. In Fig. 2, the rv. n 
program is dimensioned for up to 150 nodes, 90 points, 100 s oner: s, 
and a compressed matrix as large as 30 by 30. If tr.e wire ten: . 
makes no contact with the ground plane, the compressed mat?- x wi .' 
be exactly half as large as the full matrix. Otherwise NCM is so * - 
wr.at larger than N/2. 


c INClunr ANTI jCkiTUT : FXPJ ; 1DANT; IFFlD; I SORT; ! W i kF : />F ; ZGNM; 7GS OOl 

C INCl.lMH PD1SS;7SURF ^ 00? 

f. THlM-wlkf fiMcNNA 1 V 1 •< PFM Ff.T GROUND P L A N t 003 

C S I NUSO I UAL -f'Al r RK IN Fk » OUiNf V-OUMA I N 004 

C PROGRAM n«ir.INA?M) »«V J. rt. RICHMOND* OH I Ll S T A 1 f- UNIVERSITY 005 

C1IKIMJX IU.Yll*7ll f /H 000 

CU^Pi U C ( 30* U ) » oov 

CCIMPLPX C.l ( ISO i ( 1 SO ), U( 1 NO) , FPP( ISO) , FT 1 ( ISO) . V J ( 1 SOI 008 

COMPLFX V '» ( I SO ) ,7 l Ml SO) 000 

DIM F NS 1 I'M X 0 | 'Hi ) , YC I Of) ) , 7 C ( 90 ) , X ( wO ) • Y ( SO) ,7(90) 010 

DIMENSION IH 10 5) , UC( 100) oil 

DIMF-NMON 1 A ( liMU, I »U 100) ,'*.•»( ) no, 4 ) « NO ( 100) , C DK ( ICO) ,SDK( 100) 01? 

DIMFNN I nr., 1M 1 V> ) , I ? U SO ) , 1 3< 1 SO ) , J A ( 1 SO ) , JM ( 1 SO ) 013 

DATA P I , T P / 3 . 1 M S9 . 6 . ? 8 3 1 8 / 014 

1 F ( ?KA U HV , • JPp r • . 1 S ♦ SX , • MA X = • , 1 s , 5 X , • M I N = • t 1 S , 5X , 1 N = • , 1 S , SX , 0 1 S 

? ■ NC l* ~ 1 . 1 S ) 016 

? F I JKM A 1 t * X , • AL - • ♦ F V ! . 6, S X , 1 C MM - ■ , F H • A , SX , • FMC - * « Ml • ? ) 0 1 7 

3 FORM AT ( «*. « f- FT * . f r. ?, 3X * ' V | l e « , ?f- 1 0. 6 , 3X f ‘ 7 J I = • , ?F8. ? ) 01 R 

«♦ FORMA 1 ( 8* , • Pm • , * t , » TH » , rtx, • OOP » , 7 X, » l)B T ' , 7 X , 'GPP • , /X « » GT 1 • ) 019 

5 FORMAT ( 1 -«o ) 0?0 

6 FORMAT ( 1 X , ?F 10. 0 , A I- \ 0. ? ) 0? 1 

1CC-30 0?? 

ICJM60 0?3 

INPr^O 0?4 

INS-loo 0?S 

C THF GHlMtlKY t It- P-r ThIN-WIKF STRUCT UK F IS SPFC1FJFD IN SU?5 . IWI^F 0?6 

call I W J R F { 1 A , I K , ! NP , I NS , I N 1 , I WKC J , 1 v'K I T E , NL D, NP , NS » NkP , O? 7 

?NKS*NPGP,NSGi\A1 ,C» ? , UPH.FkC* sCAl h , T h* VG , XL * YC . /C * ZLD) 0?H 

1 1 { NSGP. IT. MPGP ) ll TfJ 5 >0 0?9 

lFtNS.GT.lNS .Ilk* NP.GT. INPJGri ll) 500 030 

r a i | ! 5.0^ ■»!! ^ ,! R# Iff. % ! CJ ,! N'S. Ik'R ! Ti * I 1 • I ?• I 3* JA * JP.* 031 

?MAX * M 1 N. M.,) t N, \C- M « N'l, NP , NS . NR P ♦ NR S . NPGP * N SOP * DC » XC » YC • 7 C ) 03? 

c NCM - SIZF OF C.MMVRtSSSn MATRIX C(1*J) 033 

JPPrNC.M-NRGP 036 

WK IT F t 6, I ) JPP, MA X, Ml N, N f NCM 03S 

WR 1 Tf t 6. 5 ) 036 

AK = TP<Al 037 

WAVH=300*/P«C 038 

WR1 H ( 6. 2 ) Al. • CM** F C 039 

WR IT H /.,*») OAO 

T P l - SC A L f T T* / a A V M O'* 1 

irtN.Lfc.O .OR. N . G I . I C J ) G U TU SOO 04? 

IFI NCM. GT. ICC )GU l»i SOO 04 3 

IFtMAX.LF.O .liK. MlN.lt.OGO Til SOO 044 

DO 9() J = l, MS 04 S 

90 {)( J >=TPL*OCt .1 ) 046 

DO 100 I-ltNP 047 

X( 1 ) = TPL*XC< 1 ) 048 

Y ( I ) -1 PL *VC ( l ) 049 

100 7(1 )=TPL-ZCt I ) OSO 

CAl L 1 DAM ( I A* 18, ICC, INS, I N T , II, I ? , 1 3 , J A , Jrt , JPP , M(), N , NCM* 001 

2ND, M D, NP, NPGP , N < S . NS * AK * C* CMM, 0, FMC , CDK , SDK* X « Y ♦ 7 . Zm, Z L D) OS? 

IF { N. F 0. 0 ) GO 10 SOD OS 3 

Il?=l 0S4 

CAU AN I I ( I A, I 8, 1 1 , I ?, I 3, 1 WKCJ. ISR II F , I 1 ?, ICC, INS, JA, JM, OSS 

? JPP, f.D, N.NC M, M ) « N t. N.NPGP * MRS, NS.C, CDK, SDK *C J * C.MM, U, LFF , G , VG, VJ , 0S6 

3Y 11 , / 11 * 7H, /LD ) 0S7 

1 F ( 1 1 ?.NF. 12 )GlT TO SOO OSH 

WK l T i ( 6, 3 ltht , Y [ I , / I I 059 

IFtG.iD.l. .AND. fc t h . FC). 0. )(»0 TU SOO 060 

WP I 1 I: t 6 * S ) 061 

L 1 MM. 5*360. /IJPH 062 


Fig. 2. The MAIN computer program. 
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WR I 7 F ( 6, 4 ) 

00 300 1 1'H: l . L 1 M 
P»«a|)PHc ( 1 f'H~ 1 ) 

f AL L If- f L IMI A t 1 8 « I NS » 1 1 t I ? * I 3 , MU, N » NO , N'S , CtK » C J * 0* 

?i pp,m , * ph, f r h, c.opp.gtt, ph. son. in, x, y, z ) 

DflP-.O 

DBT-.O 

1 f- ( GPP . r. 1 . n. ) (mp- 1 c. * A I 00 1 (I I GPP ) 

If ( C.T T. f,T . n. J »#K 1 = 10. V A LOG I CM 6TT ) 

300 Wk n |( 6, 6 I pH, I 0 , MHP . OH I , f»PP. GT T 
600 CAl L F.X I T 
F NO 


063 

064 
066 
C66 

06 7 
068 
06V 

07 0 
07 1 
07/ 
0"' 3 
0 74 


Fig. 2. The MAIN computer program - continued 
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X(I),Y(I) and Z(I) denote kx, ky and kz for point I, where 
k = If calculations are desired for a given antenna at several 

frequencies, the frequency DO LOOP will begin just below the call to 
ISORT. 


The main program calls subroutine IDANT to generate the com- 
pressed open-circuit impedance matrix C(I,J). Then subroutine ANTI 
is called to obtain the current distribution CJ(I) and the radiation 
efficiency EFT. ANTI also calculates the complex power input to the 
antenna, denoted by Y 1 1 , and the time-average input power G. If the 
antenna has only one generator and VG(JGN) = (l.,0.), then Y 1 1 and 
Zll denote the antenna admittance and impedance, respectively. 

Finally, the antenna pattern is obtained by calling subroutine 
IFFLD. TH and PH denote the spherical coordinates c and $ (in 
degrees) of the distant observer. GPP and GTT denote the {-polarized 
and ''-polarized power gains, respectively, and DBP and DBT are the 
decibel versions. The user may want to increment o as well as < , 
but this will require only a trivial change in the main program. 

IFFLD is called once for each look angle (o,0- When o = 90°, GPP 
will vanish if the program has set up a valid system of images. . 


IV. AN EXAMPLE 


Fig. 3 shows a simple antenna and its image, with a dotted 
line to indicate the ground plane. In Fig. 1, subroutine IWIRE sets 
up the following input data for this antenna: 


VG(4) 

= 1 

NLD = 

0 

CMM = 

1 . 

NSGP 

= 2 

NPGP 

= 2 

NRS = 

4 

NRP = 

5 

is planar 

antenna 


unit voltage generator at endpoint IA of segment 4 
no lumped loads 

the wire conductivity is 1 megamho/m 
2 real segments connect to the ground plane 
2 points on the ground plane 

4 real segments 

5 real points 


and the numbering system is shown in Fiq. 3. 


8 segments (NS = 8), 
Note that the lowest 


numbers are assigned to the two points on the ground plane and the 
two segments terminating on the ground plane. The points and seg- 
ments must be labeled with consecutive positive integers 1,2, 3, . 
For a given segment J, it makes no difference which end is labeled 
I A ( J ) . In Fig. 3, each numeral located near a dot is the index I 
of that point. Each numeral located near the center of a line is 
the index J of that segment. 
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Fig. 3. Points and segments on a simple wire antenna 
lying in the xz plane. 


Fig. 4 shows the same antenna ana the eight dipole modes 
defined by subroutine ISORT. The arrows indicate the reference 
directions for the mode currents and voltages. The mode index 
number I is placed near the terminal point 12(1). Mode I is a 
sinusoidal basis function which vanishes at the endpoints 11(1) 
and 13(1) and has unit current at the terminal point. These are 
overlapping subsectional bases, and mode I extends over two inter- 
secting segments JA(I) and JB ( I ) . The reference direction for node 
currents and voltages is from II to 12 toward 13. In Fig. 4, 
inodes 1 and 2 have terminols at the ground plane, with segment JA 
above and segment J3 below the ground plane. This type of node has 
no image. Modes 3, 4 and 5 have images. The size of the compressed 
matrix is NCM = 5. If we did not take advantage of the ground-plane 
symmetry, the matrix size would be N = 8. 

Table I presents some of the output data for this example 
and Table II lists the elements of the compressed matrix C(I,J) on 
return from subroutine IDAfiT. From Table I, the calculated imped- 
ance is Z 1 1 = 959 + j 664 ohms. For the same antenna with perfect 
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4 



Fig. 4. Mode map for the antenna shown in Fig. 3. 


conductivity, Zll = 879 + j 749 ohms. The calculated results 
should not be considered accurate without checking the behavior as 
the wires are subdivided into more segments. The longest segment 
should not exceed A/4 . The thin-wire and delta-gap formulations 
are justified most readily if the wire radius does not exceed 0.007A. 
Fortuitously, satisfactory results are often obtained for closed 
wire loops even when the wire radius is as large as 0.02A. For 
dipoles, an upper limit of 0.007A is recommended. 

With 20-ohm resistive loads inserted in each end of each 
segment of the antenna in this example, the calculated impedance 
is 710 + j 206 ohms and the efficiency is 45.5 per cent with CMM = 1. 
With and without lumped loads, the results obtained with the present 
computer program show satisfactory agreement with those obtained 
with the program in reference 1. (A new version of subroutine SCANT 
was used for tnese tests because the original version in reference 1 
does not always handle lumped loads properly.) 
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TABLE I 


1NP 

=> 93 

NP = 

e 

INS 

= 1 00 

NS = 

C 

J 

IA< J) 

1 B ( J > 



1 A ( X ) 

i n c k ) 

DC < J ) 

1 

1 

3 


5 

i 

6 

. s 0 :• r o 

2 

2 

4 


6 

2 

7 

.50000 

3 

3 

4 


7 

6 

7 

1 .00000 

4 

3 

s 


6 

6 

6 

• 99985 


1 

( I > 

c I ) 

2CC i ) 

J 

XC ( J ) 

YC C J > 

ZC < J > 

1 

1 ■ 

*0.00000 

0 as 





2 

0 . 00003 

0.03000 

0.3.' 030 





3 

1 .00000 

0 . 00000 

.50030 

6 

l .00000 

0.00000 

- .53000 

4 

0 .O0C93 

0 .00000 

.50000 

7 

0.03330 

0 .03003 

- .533 0 3 

5 

1 .73730 

0.00000 

1 .20730 

8 

I .70700 

0 .0000 0 

- 1 .20720 

I 

JA JB 

I 1 12 

I 3 


K J A 

JP I 1 

12 13 


1 

: 

r t 

3 

1 

6 







2 

2 

6 

4 

2 

7 







3 

1 

3 

1 

3 

4 

6 

5 

7 

1 

6 

7 

4 

3 

4 

4 

3 

5 

7 

7 

8 

7 

6 

8 

5 

2 

3 

2 

4 

3 

3 

6 

7 

2 

7 

6 


JPP* 

3 

MAX a 

3 MIN = 

1 N* 

8 

NCM = 

5 

AL* 

.0001 00 

CMM* 

1 .0330 

FMC « 75.00 




J= 

4 

VG ( J ) * 

1 .00 

0.00 





1 

MAGNITUDE. 

PHAGE 

PEAL 


I MAG I r.’APY 

2 

3 

4 

5 

EFF «= 
Pil 

1 .033 
.873 
.730 
.271 
.638 

90.44 Y1 1 

TH 

99.9 -.0035441 

91.1 -.0003597 

-77.6 .0034892 

-34.7 .0037048 

-85.1 .0331362 

= .033705 -.000488 

D3P D3T 

Z1 1 

GPP 

.03311 £3 
.0027537 
- . 0 2 * 2 r » 7 o 

-.P004eso 

-.0023153 

= 959.07 664.07 

fVT 

O • 

85- 

0.03 

.66 

0.00 

1.16 

20. 

85. 

-47.44 

1 .24 

.00 

1 .33 

40. 

85. 

-39.49 

2.63 

.00 

1 .83 

60. 

65. 

-33.35 

4.17 

.00 

2.61 

80. 

85. 

-29.46 

5.31 

.03 

3.40 

1 00 . 

65. 

-27.93 

5 .77 

.03 

3.77 
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TABLE II 


Compressed Impedance Matrix 


I 

J 

C(I,J) 

C(J, 

*) 

1 

1 

16.1 

"j 

720.5 

1671 

-j 

720.5 

1 

2 

7.6 

“J 

6.7 

7.6 

-j 

6.7 

1 

3 

-15.9 

-j 

1059.9 

- 7.9 

-j 

529.9 

1 

4 

-18.3 

“J 

16.8 

- 9.4 

-j 

8.4 

1 

5 

- 5.2 

+ j 

83.0 

- 2.6 

+j 

41.5 

2 

2 

16.1 

“j 

720.5 

16.1 

-j 

720.5 

2 

3 

- 5.2 

+ J 

83.0 

- 2.6 

+j 

41.5 

2 

4 

-12.7 

“J 

77.3 

- 6.4 

-j 

38.6 

2 

5 

-15.9 

-j 

1059.9 

- 7.9 

-j 

529.9 

3 

3 

21.2 

-j 

326.0 

21.2 

-j 

326.0 

3 

4 

- 9.2 

"J 

22.6 

- 9.2 

-j 

22.5 

3 

c 

- 9.2 

“j 

396.0 

- 9.2 

-j 

396.0 

4 

4 

51.3 

+ j 

60.8 

51.3 

+j 

60.8 

4 

5 

22.1 

+ j 

420.9 

22.1 

+j 

420.9 

5 

5 

21.2 

"j 

326.0 

21.2 

-j 

326.0 


V. SUMMARY AND CONCLUSIONS 

A computer program is presented for a thin-wire antenna over 
a perfectly conducting ground plane of infinite extent. The analysis 
is performed in the frequency domain, and the exterior medium is free 
space. The antenna may have finite conductivity and lumped loads. 

The output data includes vhe current distribution, impedance, radia- 
tion efficiency and gain. The program uses sinusoidal bases and 
Galerkin's nethod and takes advantage of the ground-plane symmetry 
to reduce the storage requirements and computation costs. The 
subroutines are included in alphabetical order in the Appendices with 
a brief expl anation . 
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APPENDIX 1. Subroutine ANTI 


Subroutine ANTI is listed in Fig. 5. Between statements 14 
and 30, this subroutine sets up the excitation voltages CJ ( I ) and 
VJ(I) with the aid of the delta-gap model and the input data for the 
generator voltages VG(J). ANTI calls CROUT to obtain a solution for 
the simultaneous linear equations . On return from CROUT, the dipole 
mode currents are stored in CJ(I). The image currents are stored 
in CJ(K) in the DO LOOP ending with statement 80. The DO LOOP 
ending with statement 90 calculates the complex power input Yll and 
the time-average power input G. The power dissipated (DISS) in the 
lumped loads and the imperfectly conducting wire is obtained by 
calling PDISS. Finally, the radiation efficiency EFF is calculated. 

If IWRCJ is positive, ANTI writes a list of the dipole mode 
currents CJ ( I ) . This list includes the normalized current magnitude, 
the phase in degrees, and the real and imaginary parts of the 
current. 


APPENDIX 2. Subroutine CROUT 

CROUT, listed in Fig. 6, solves a system of simultaneous linear 
equations with complex coefficients. This subroutine uses the method 
of P. D. Crout. Although this subroutine does not use pivoting, it 
is efficient and accurate in the present application. The input 
data are defined as follows: 

C ( I , J ) complex coefficients in the simultaneous equations 

S(I) excitation coluim 

ICC dimensions of C and S 

ISYM zero or one for symmetric o v ' nonsymmetric matrix 

IWR one or zero to write or suppress the solution 

112 one or two if C is original or auxiliary matrix 

N size of the square matrix C 

Of course, N must not exceed ICC. If IWR is a positive integer, the 
solution will be printed out with the following definitions: 

I index number of the solution S(I) 

SNOR normalized magnitude of S(I) 

SA absolute magnitude of S(I) 

PH phase of S(I) in degrees 

On the first call to CROUT, C ( I , J ) contains the original 
matrix. 112 = 1 and CROUT generates the auxiliary square matrix, 
overlaying it in the sane location C and destroying the original 
matrix, then CROUT proceeds to generate the solution, storing it in 
S(I) and destroying the original excitation coium. 
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Next we might want another solution of th ame system of 
simultaneous linear equations but with a new ex- ition colurm. This 
could be obtained by recalculating the original matrix C(I,J) and the 
new excitation column and 'calling GROUT again with 112 = 1. However, 
there is no need to recalculate C(I,J). Instead generate the new 
excitation colunn, set 112 = 2 (or any integer other than 1) and call 
CROUT again. CROUT uses less computer time when 112 differs from 1. 


APPENDIX 3. Subroutine EXPJ 


Subroutine EXPJ, listed in Fig. 7, evaluates the exponential 
integral defined as follows: 


W12 


V2 

e - V - - Y - = Ej (VI) - Ej(V2) + j 2n* 
v 

VI 


where the integration path is the straight line from VI to V2 on the 
complex v plane and 



E j ( z / denotes the principal branch of the exponential integral. To 
generate W12, EXPJ calculates Fj(Vl), subtracts E 1 ( V2 ) and adds j2nu. 
The integer n is zero unless the straight-line integration [ .ch 
intersects the negative real v axis at a point between VI and V2. 

When there is such an intersection, n = 1 if VI lies in quadrant 1 or 
2 and n = -1 if VI lies in quadrant 3 or 4. 
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APPENDIX 4. Subroutine IDANT 


Subroutine IDANT is listed in Fig. 8. This subroutine stores 
the quantities CDK(J) = cos kd-j and SDK(J) = sin kd.; where dj is the 
length of segment j. The program wri tes AK, DMAX and DMIN and aborts 
if 


a. the length of the shortest segment is less than the wire 
radius, or 

b. the longest segment has a length d such that kd exceeds 3, 
or 

c. the wire radius a is such that ka exceeds 0.1. 

IDANT calculates the elements in the compressed impedance 
matrix C(I,J) as follows. Select a source segment K and a receiving 
segment L, where K and L range from 1 to NS. The mutual impedances 
Pll , P12, P21 and P22 between the two segments are obtained by calling 
ZGMM if K = L, ZGMM if the segments intersect, or ZGS if segments K 
and L do not intersect. 

Nav select a test dipole I sharing segment K, and an expansion 
dipole J sharing segment L. Add the appropriate segment-to-segment 
impedance to the dipole-to-dipole impedance C(I,J). When this 
procedure has been completed at statement 200, the impedances C(I,J) 
are appropriate for a perfectly conducting thin-wi re system with no 
lumped loads. 

Between statements 200 and 262, the impedance matrix C ( I , J ) is 
modified to account for the finite conductivity of the wire antenna. 

The surface impedance ZS is obtained by calling ZSURF. For each 
segment K, the program selects a test dipole I and an expansion dipole 
J sharing this segment. The contribution to C(I,J) associated with 
finite conductivity is ZSAM if dipoles I and J have terminals at the 
same end of segment K, and Z0PP if they have terminals at opposite 
ends. C(I,J) is not affected unless dipoles I and J share one or 
two segments. 

Between statements 262 and 280, the irpedance matrix is modified 
to account for the lumped loads. Each diagonal element C(I,I) is 
modified by adding the impedance of the lumped load inserted at the 
terminals of node I. If modes I and J share a segment and have termi- 
nals at the same point, C(I,J) is modified by adding or subtracting 
the impedance of the lumped load inserted at the terminal end of 
this segment. (Add or subtract ZLL) if mode currents I and J have the 
same or opposite reference directions on the shared segment.) 
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APPENDIX 5. Subroutine IFFLD 

Subroutine iFFLD, listed in Fig. 9, calculates tht far-zone 
field of the thin -wire antenna. 

Let (r,P,0 denote the spherical coordinates of the distant 
observer, ano let E 0 (I) and F.(I) denote the electric field intensi- 
ties of dipole inode I with unit terminal current. Then 

EPP(I) = (r/x) e jkr E ^ ( I ) 

ETT(I) = (r/x) eJ' kr Eg ( I ) 

The field of sinusoidal dipole mode I may be regarded as the sum of 
the fields of each of its two segments. The field of segment K is 
obtained by calling subroutine 7.FF, and EPP(I) and ETT(I) «re 
generated by adding the appropriate numbers obtained from two 
different calls to ZFF. In the DO LOOP ending with statement 260, 
the antenna field is calculated as a weighted sum of the mode fields 


CJ(I) EPP(I) 


C0(I) ETT(I) 


where CJ(I) denotes the terminal current of mode I and EPH and ETH 
denote the dimensionless range-independent form of the antenna fields 
E 4 and E 0 . 

G denotes the time-average input power to the antenna, and 
GPP and GTT are the ^-polarized and O-polarized power gains. Sub- 
routine IFFLD is called once for each angular direction. In the 
input data supplied to this subroutine, PH and TH denote and e in 
degrees . 

This subroutine is useful for wire antennas with or without a 
ground plane. For an antenna over a ground plane, IFFLP must be 
supplied with information on the couple te system including the image. 


as follows: 

N 

EPH = l 
1 


N 


ETH = 
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APPENDIX 6. Subroutine I SORT 


Subroutine ISORT, listed in Fig. 10, is described briefly in 
Section III. This subroutine sets up the image segments and points 
and calculates the segment lengths. Then it checks the input data 
for consistence. The data are considered inconsistent and the run 
is aborted i f 

a. NPGP is greater than zero but one of the points on the 
ground plane has no segment (with index J l^ss than or 
equal to NSGP) connected to it, or 

b. a real point situated above the ground plane has no 
segment (with index J less than or equal to NRS) 
connected to it. 

Between statements 32 and 50, this program calculates the 
number of es N on the complete structure. The run is aborted if 
the dimensions are inadequate. 

Between statements 50 and 58, the program sets up the nodes 
that v.ill not have images. The number of node s of this type is 
NPiiP* and these modes have the lowest index numbers. Mode I has 
terminal point 12(1) = I on the ground plane, endpoint 11(1) is above 
the ground plane* endpoint 13(1) is the corresponding image point 
below the ground plane, and segment JA(I) is the lowest -numbered 
real segment with endpoint I. 

Between statements 58 and 65, the program sets up the rest of 
the real modes. Modes of this type have the terminal point 12 on 
or above the ground plane. Each of these real modes (with index I 
qreater than NPGP) has an image which is established between 
statements 65 and 75. 

Below statement 75, the last part of the program counts the 
number of dipole modes sharing segment J, denoted by ND(J). It also 
stores a list of the dipole modes sharing segment J, denoted by 
MD(J,K). A segment may be shared by as many as four modes. 


APPENDIX 7. Subroutine PDISS 

Subroutine PDISS is listed in Fig. 11. This subroutine 
calculates the time-average power (DISS) dissipated in the lumped 
loads and the imperfectly conducting wire. The power is calculated 
for one segment at a tiro, and the total power dissipated is the sum 
of the powers dissipated on the various segments. On segment K, CJA 
and CJB denote the currents at endpoints I A ( K ) and IB(K). RLA and 
RLB denote the lumped resistors inserted in segment K at endpoints 
IA 3nd IB. 
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This subroutine is suitable for a wire structure in free space 
and also for a wire structure over a perfect ground plane. If there 
is no ground plane, the total number of segments NS mu^t be supplied 
as the tenth calling parameter, and DISS denotes the power dissipated 
on the entire structure. If there is a ground plane, the number of 
real segments NRS is supplied instead, and DISS denotes the po wer 
dissipated on the real segments. 


APPENDIX 8. Subroutine 7.FF 

Subroutine ZFF, listed in Fig. ;2, calculates the far-zone 
field of a sinusoidal electric monopole in free space. Tne monopole 
has endpoints at (XA,YA,ZA) and (XB,YB,Z8). (These symbols denote 
kx, ky and kz.) Let E denote tne electric field intensity. The 
di mensiort less range-independent field is defined by 

F = (r/\) e-i kr E 

FP1 and FT1 denote F. and Jr r for the mode with unit current at 
(XA,YA,ZA). EP2 and’ET2 denote F, and F, for the node with unit 
current at (XB,YB,78). The far field vanishes in the endfire dire 
where GK - 0. 


APPENDIX 9. Subroutine ZGMM 

Subroutine ZGMM, listed in Fig. 13, calculates t he .mutual 
impedance between two filamentary monopoles with sinusoidal < , 
distributions. The dipole-dipole mutual impedance in Fg. 20 r 
reference 2 is the sum of fc* r monopole -monopole mutual impede 
The nonopole impedances are calculated by ZGS with Simpson's r 
by ZGMM with eloseo-ferm expressions in terms of exponential . 

If the monopoies are parallel , let the z axis be paralle 
both ronopoles. The coordinate origin nay be selected arbitral 
SI and S2 denote the z coordinates of the endpoints of the te5 * 
nonopole, T1 and T2 are the z coordinates of the endpoints of t 
expansion monopole, and D is the perpendicular distance (displat 
between the monopoles. The mutual impedance of parallel r onopo. 
calculated in the last part of ZGMM below statement 110. 

For skew monopoles, let the test monopole s lie in the xy 
plane and the expansion noncoole t in the Plane i = D. (D is the 
perpendi cular distance between tne parallel planes.) If the rone: 
are viewed along a line of sight parallel with the z axis, tne ex- 
tended axes of tne two ronopoles will appear to intersect at a pon 
on the xy plane. Let s measure the distance along tne axis of the 
test monopole with origin at the apparent intersection. SI and S2 


denote the s coordinates of the endpoints of the test monopole. 
Similarly, let t measure distance along the axis of the expansion 
monopole with origin at the apparent intersection. T1 and T2 denote 
the t coordinates of the endpoints of the expar.s in monopole. Let 
s and t be unit vectors parallel with the positive s and t axes, 
respectively. Then CPSI = s • t = cos \^. The monopole lengths are 
d s and dt, and the remaining input data are defined as follows: 

CGDS = cos kd s 

SGD1 = sin kd s 

SGD2 = sin kdt 

ZGMM calls EXPJ for the exponential integrals. ZGMM is 
specialized for sinusoidal monopoles in free space. In ZGMM the 
input data SI, S2, T1 , T2 and D denote ks , , ks., , ktj , kt 2 and kd, 
respectively. Otherwise, ZGMM is the same as GGMM. 

The output data from ZGMM are the impedances PIT, P12, P21 
and P22. In defining these impedances, the reference direction is 
from SI to S2 for the current on monopole s, and from T1 to T2 for 
the current on monopole t. In the impedance P-h , the first subscript 
is 1 or 2 if the test dipole has terminals at SI or S2 on monopole s. 
The second subscript is 1 or 2 if the expansion dipole has terminals 
at T1 or T2 on monopole t. The endpoint coordinates SI, S2, T1 and 
T2 may be positive or negative. The monopole lengths d s and dt are 
assumed positive in defining the input data CGDS, SGD1 and SGD2. 

For parallel monopoles, CPSI = 1 or -1. SI, S2, T1 and T2 
are cartesian coordinates for parallel monopoles and spherical 
coordinates for skew monopoles. For skew monopoles, the radial 
coordinates SI, S2, T1 and T2 tend to infinity as the angle ^ tends 
to zero or -n . Therefore, if the monopoles are within 4.S degrees 
of being parallel, they are approximated by parallel dipoles. 


APPENDIX 10. Subroutine ZGS 

Subroutine ZGS, listed in Fig. 14, calculates the mutual 
impedance between two filauentary monopoles with sinusoidal current 
distributions. (The dipole-dipole mutual impedance in Eq. 20 of 
reference 2 is the sum of four monopole -monopole mutual impedances.) 
The endpoints of the axial test monopole s are at (XA,YA,ZA) and 
(XB,YB,ZB), and the endpoints of the expansion monopole t are at 
(X1,Y1,Z1) and (X2,Y2,Z2). DS and DT denote the lengths of mono- 
poles s and t. Dimensionless forms are used for the input data. 

For example, XA, AK, DS and DT denote kx a , ka, kd s and kdt- CAS, 

CBS and CGS are the direction cosines of monopole s, and CA, CB and 
CG are the direction cosines of monopole t. 
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If INT = 0, Z6S calls ZGMM for the closed-form impedance 
calculations. Otherwise ZGS calculates the mutual impedance via 
Simpson 1 s-rule integration with the following number of sample 
points: IP = INT +1. If the nonopoles are' parallel with small 

displacement, ZGS calls ZGMM to avoid the difficulties of numerical 
integration . 

For' the fields of the test monopole, ZGS uses Fqs. 75 and 76 
of reference 2. The current distribution on the expansion monopole 
is given by Eq. 74 of reference 2. With an origin at (X1,Y1,Z1), 
the coordinate T measures distance along the expansion monopole. 

Thus T is the integration variable. 

Let the coordinate s measure distance along the test monopole 
with origin at (XA,YA,ZA). From any point T on monopole t, construct 
a line to the test monopole such that the line is perpendicular to 
the test monopole. SZ denotes the s coordinate of the intersection 
of this line with the test monopole. The length of the line is the 
radial coordinate p, and RS denotes . 2 . R1 and R2 are the distances 
from (XA,YA,ZA) and (XB,YB,ZB) to the point T. Cl is the current 
at T for the mode with terminals at (X1,Y1,Z1), and C2 is the current 
at T for the other mode with terminals at (X2,Y2,Z2). C denotes 
the Simpson 's-rule weighting coefficient. 

Below statement 300, ZGS performs some analytic geometry in 
preparation for calling ZGMM. The remaining part of this Appendix 
concerns this last part of subroutine ZGS. 

Lgt s denote a unit vector from (XA,YA,ZAl toward (XB,YB,ZB), 
and let t denote a unit vector from (X1,Y1,Z1) toward (X2,Y2,Z2). 

Then s • t = cos ^ = CC where ^ is the angle formed by the axes of 
the two monopoles. Let monopole s lie in one plane P s and monopole t 
in another parallel plane Pt. CAD, CBD and CGD are the direction 
cosines of the unit vector d = t x s / sin ^ which is perpendicular 
to both planes. To obtain the distance DK between the planes, we 
construct a .vector R, , from (XA,YA,ZA) to (X1,Y1,21) and take 
DK = R n • 1 

Construct a line from (X1,Y1,Z1) to the test nonopole, such 
that the line is perpendicular to the test monopole. SZ denotes the 
s coordinate of the intersection of this line with the test monopole, 
and the cartesian coordinates of this intersection are XZ, YZ and ZZ. 
The direction cosines of ? x fl are CAP, CBP and CGP. 

From the point (XI ,Y1,Z1) in plane P^, construct a perpendicu- 
lar line to the point (XP1 , YP1 ,ZP1 ) in the plane P s . This line is 
parallel with d and has length DK. Let R^ represent a vector from 
(XZ,YZ,ZZ) to ( XP1 ,YP1 ,ZP1 ) . PI denotes R-(s x d) . SI and T1 are 
defined in Appendix 9. 
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Subroutine Z6S is essentially the sanva as 66S except the 
medium is specialized to free space in ZGS. 


APPENDIX 11. Subroutine ZSURF 

Subroutine ZSURF, listed in Fig. 15, calculates the surface 
impedance of a solid circular-cylindrical wire with exterior exci- 
tation. ZS denotes the surface impedance in ohms, and the input 
data are defined as follows: 

AK ka, where a is the wire radius 

CMM conductivity of the wire in megamhos/m 

FMC frequency, MHz 

The surface impedance is defined by Z s = Ez/Hi where the fields E an u 
II are evaluated at the surface of the wire. This subroutine calcu- 
lates the impedance for the lowest order cylindrical mode with fields 
E z and H. independent of <i . The wire is considered to be a good 
conductcV in the sense that the displacement current is negligible 
in comparison with the conduction current. In the present applica- 
tion, we require the surface impedance appropriate for the current 
distribution I(z) = sin kz. For a highly conducting wire, however, 
this impedance is considered to be the same as that for a uniformly 
distributed current. BER, BEI, BERP and BEIP denote the Kelvin 
functions ber, bei, ter' and bei ' with argument x. When x is less 
than 8, BES and BES1 denote the Bessel functions J 0 and Jj with 
argument 

z = x e’j"/ 4 • 


When x is greater than 8, 


.707x 

J n /BES = J , /BE SI = £ 

* / 5ttx’ 
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suoftnur inf ant i i i a, ih, n , 12, l iwkc j, i write* i 12, icc, ins, ja, jb, 001 

2JPP,M!),N,NCM f NO,Nl.|i,NPGP, NkS, \ S . C , CDK , SDK • C J , CMM , 0 , E F F • G , VG • V J , 002 

3Y11 , 7 It . /*, ZLD) 003 

complfx cncc . icr.i *cj( i ) , vc,i n , vju > »ZLom 004 

COMPLEX Y1 1 t 1 1 1 , Zh, rj\ . Vj! , fvv 005 

DIMENSION 1 All) , Ml if, 1 ),CUK( 1 ) . SOK (1) . Dll) , i 1 (1) , 1 2 ( 1 1 , I 3 (1) 006 

D IMF NS I U J MD( INS, 4) , JAl 1 ), Jh( 1), IGfcNI 1 ) , JGLN{ 1 } 007 

1 FORMA T ( HX, M= « , l >t 5<, 1 Vbl J ) = 1 , 2F10.2 1 00ft 

2 FORMAT ( rtX, • j = ' t I :>, sx, • Z L 0( J)- *, 2F 10. 2) 009 

3 FORMAT I 10X, •!•./*/., ' f»AGN I I UUF ' , 3* , • PHASE • ,VX, ' s F A L ' ,8X, • IMAGINARY' ) 010 

4 FORMA 1 ( IX, 1 | 10, 1 F10. 3, 1 F10. 1 , ?F 15. 7 1 Oil 

6 FORMAT (1H0) 012 

t F f II 2.GT.0) 1 12 = l 013 

1 F C ! WR IT F.LF. 0)1*0 llJ 14 01 A 

0(1 10 J = 1,\RS 016 

AVG=CAftS( VG( J )) 016 

10 I F ( AVG. Gl.0.01 ) WR 1 7 E ( 6, 1 1 J, VG( J ) * 017 

WR IT F ( 6, 6 ) 01 H 

IFINl 0.| b.O)GO 111 14 019 

7 M AX = . 0 020 

OD 12 j= 1 , NR S 021 

A7L=CARS( 7LIX .1 ) 1 022 

IF I A71. IT .0. 1 )GU 11 12 023 

2 M A X = 1 . 0 ?4 

WRITf (6,2) Jf 7UM J ) 025 

12 CONTINUE 026 

IF(2MAX.GT.0,6)WRI1E(6,6) 027 

14 00 30 I = 1 , N C M 028 

CJ( I ) = ( • 0, . 0 ) 029 

L =2 030 

FAC-1. 031 

K = J A ( 1 ) 032 

I F ( 1 . GT . N PGP ) GO TO 15 033 

L = 1 034 

F AC = 2. 035 

1 h ( JH ( I ) • L T • K IK = Jft ( I ) 036 

15 00 25 .< K - 1 , L 037 

KA= I A( K ) 03B 

K R = I ft ( K ) 039 

J J = K 040 

FI=FAC 041 

IF (KB. FO. I ?( 1 ) 100 TO 22 04? 

IF( KR. FO. ! 1 I I) )F 1 =-FAC 04 3 

cji i )=cj< i )♦: i-vgi jj) 044 

GO TO 26 045 

22 1F(KA.F0.I3( 1 ))FI=-FAC 046 

JJrK+NRS 047 

CJI 1 )=CJ( I ) ♦F 1 * VG ( JJ ) 048 

25 K=JR( 1 ) 049 

VJ( I ) -CJ( I ) 050 

Kef+JPP 051 

30 1 F ( 1 . GT . NPGP 1 V J ( K) s-V J ( I ) 052 

ISYM=1 053 

lF(N.EO.NPGP) I S YM= 0 054 

I WR = 0 055 

CALL CROUT(C,CJ, ICC, 1SYM, 1WK, I 12.NCM) 056 

112=1? 057 

CMAX=.0 058 

00 HO 1=1, NCM 059 

CA=CARS(CJ{ 1 ) 1 060 

K= I ♦ J PP 061 

1 F ( 1 .GT.NPGP)CJ(K) = -C J ( 1 1 062 


Fig. 5. Subroutine ANTI. 
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80 !F(CA.r>T.CMAX)OlAX = CA 063 

1F( 1WRCJ.GF. 1 }WR1 TM6, 3) 064 

Yll=(.0,.0l 065 

G-l. 066 

FFF=.0 067 

711 r C • 0 * *0) 068 

1F(CMAX.LF.0# I GO TO 500 069 

g~ . o • 070 

on 90 I -1 , N 071 

CJl=GJ( 1 ) 07? 

VJI = V j ( I ) 073 

DYY=CJl*CnNJGI7Jl y 074 

I F ( I • L F • NO’ ) Y \ l = Y 1 1 4 DY Y 075 

G = G 4 P F A I { 0 Y Y ) 076 

IF( 1WRCJ. LF.OIGO T( i 90 077 

I F ( | . GT • NCM ) GO TO 90 078 

CA=CA8S(CJI ) /CrtAX 079 

PHr-.O • 080 

1MCA.GT. l.E-30>PH=57.2967H*ATAN?( A1MAG(CJI l.kEAKCJl)) 081 

wk I IF (ft, 4) l,CA t PH f CJI 08? 

90 CONTINUE 083 

III 1WPCJ.GE. 1 )WP|TL*(6,5) 084 

G=G/2. 085 

711=1# /Y 11 086 

FFFrlOO. 087 

1F(CMM.LE#0# .AMO. NL D. Lfc * 0 ) GO TU 500 088 

CALL P I) I S S ( I A t IB, INS* I 1 « I ? * I 3 * WD* NOi MLO» NKS t C J i CMM f Di CDK * 089 

2S()K,01 SStZH, 2L0) 090 

FFF = 100#*(G-MSS)/G 091 

500. RF TURN 09? 

FNO 093 

Fig. 5. Subroutine ANTI - continued 
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SUBROUTINE CRCJUKC. St ICC, ISYM, 1WK, I 12,N) 001 

COMPLEX Cl ICC, ICC).S(I) 002 

COMPLFX F,P,S$,1 003 

? FORMAT ( IX, 1 13, If 10. 3, IMS, 7, IF 10. 0) 00* 

5 FORMAT (IhO) 005 

1 F ( I 12. NF. 1 ) GO TO 22 006 

IF(N.FO. I )SI 1 ) = S ( 1 l/CI 1, 1 ) 007 

1 F ( N. FQ. 1 ) GO TO 1 00 008 

I F ( ISYM.NE.01G0 TO 8 009 

DO 6 UliN 010 

DO 6 J=I,N Oil 

6 c(j,n=cn,j) 012 

8 F =C (1,1) 013 

00 10 1=2, N 016 

10 Cll,LKll,l)/f 016 

00 20 1=2, N 016 

LLL=L~1 017 

Dll 20 1=1, N 018 

F=C(1,U 019 

DO 11 K = 1 , L L L 020 

11 F = F-C|],K)*CU,U 021 

C ( I , L ) =F 022 

IFIL.FO. 1 ) GO TO 20 023 

P=C ( L , L ) 026 

1 F ( ISYM.E0.01 GO TO 15 026 

F=C( L, I > 026 

00 12 K = 1 , L L L 02 7 

12 F = F— C ( l , K )*C(K f I 1 028 

C(L, 1 )=F/P 029 

GO TO 20 020 

16 F =C ( I , L 1 031 

C(l» 1 )=F/P 032 

20 CONTINUE 033 

22 DO 30 1=1, N 036 

P=C(L,L1 035 

T=S(L) 036 

IF1L.E0. 1 ) GO TO 30 037 

L l L = l. - 1 038 

DO 25 :<=1,LLL 039 

25 T - T -C ( L , K ) * S ( K ) 060 

30 S ( L 1 = T / P 061 

DO 38 L =2 , N 062 

]=N-L+1 063 

11=1+1 066 

T = S ( 1 1 065 

DO 35 K = I 1 , N 066 

36 T = T — C ( !,K)*S(K) 067 

38 S ( I ) = T 068 

IFMWR.LF.01 go to 100 069 

WR 1 T E ( 6 , 5 ) 060 

CM0R=.0 061 

DO 6 0 I = 1 » N 062 

S A =C AHS ( S( 1 1 ) 063 

60 IF ( S A . G T • CNI IR ) C NOR = S A 066 

1 F ( C NOR .IF .0. )CN0R = 1. 066 

DO 66 I=1,N 066 

SS=S(I) 067 

S A = C A B S ( SSI 068 

SNOR = SA/CN()R 069 

PM=.0 06 0 

IF( SA.GT .0. ) P H = 5 7 . 29578* AT AN2( A I MAGI SS), REAL (SSI ) 061 

64 WR1 TE ( 6, 2 M , SNOR, SA, PH 062 

WR I T E ( 6 , 6 ) 063 

100 RETURN 066 

ENO 063 


Fig. 6. Subroutine CROUT. 
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* SUBROUTINE FXPJ< VI, V2, W| 2 > OOl 

COMPLEX EC,E19,S. T.UC, VC, VI, V2, W12,Z 002 

0 I MENS I ON VI21 ),W(21 ),l)l 16),t < 16) 003 

DATA V/ 0.22284667E OC, 004 

20.11RR9321E 0 1 , 0. 2992 7 36 36 01 , 0. 9775) 436F 0 1 , 0. 9»374674£ 01, 009 

?0.1998?874E 02, 0. 9 330?H 1 2 K-0 1 , 0. 4V26V 1 / 4P 00,0. I 2 I 999946 01, 006 

20.22699499E 01,0. 3667*62271: 0 1 , 0. 5 42533666 01 , O. /969V 1 62E 01, 007 

20. 1012022RE 02 , O. 1 3 1 30282F 02,0.166944086 (>2,0.207 / 6 4 7 9 6 02, 008 

20.296238946 02 , 0 . 3 1 4 079 1 Vr 02,0.385306836 02 , 0. 48026086 t 02/ GOV 

DATA w/ 0.45896460E 00, 010 

20. 9 1 7O00*3t- 00,0.1)3373386 00, 0. 1 0399 1 V7E-0 1 , 0. 261 0 1 7206-03, Oil 

20.89854791F-06, 0.2 1 8 ? 3 4 8 7 E 00, 0. 34?? 101 76 00 , 0. 26 30?76>‘F 00, 012 

20. 126425826 00,0. 4020686SF-0 1 , 0. 856387 /86-0? , 0. 1212'*361E- 2, Oi3 

20. 1 1 1 6 74 4 0 6 -0 3 , 0 • 6 4 59 9 ?6 7 F - 0 5 , 0 . 2 2 2 6 3 1 69 6 - 06 , 0 • 4 ? 2 7 4 304 -08 , 0 1 A 

20. 392189736-10,0. i '.965 I 52C - l 2, 0. 148302/06-1 9,0. 160059496-19/ 019 

DATA ()/ 0.224948426 02, 016 

2 0.744U568F 02, -0. 4 1 '*3 1 9 766 03 , - 0. 7H 7 96 3398 02, 0. 11254?a<,E 02, 017 

2 0.160217616 03 , -0 . 23862 1 991* 03, -0. 9009668 /E 03 . -O. 68*878946 02, 0 1 r. 

2 0. l??S477Ht 02, -0. 10161976b 02,-0. 4?2 l 99916 01. 0.7V729681E 01, 019 

2-0.210699746 02, 0.220464906 01, 0.697262446 01/ 020 

DA 1 A F/ 0.21 1 0 3 i 0 7 E 02, 021 

2-0. 379997876 03,-0.974892206 02. . 1 29006726 03, 0.179492266 02, 022 

2-0. 1 29 i 0931 6 03,-0. 997099746 0*., . 13524801F 02, 0. 146967216 03, 023 

2 0.1794V52HE 02 , -0. 3298 1 0 1 4 1 00, . 310288366 02, 0.M69V697E 01, 024 

2 0. 222 3696 1 E 02, 0. 391248926 02, 0.81636799E 01/ 029 

2 = V 1 026 

DU 100 JlM*l f 2 027 

X = R F A L ( Z ) 028 

Y« A I MAG ( Z ) 029 

E ] 9 r ( . Of . 0 J 030 

AtUCAKS(Z) 031 

1 F ( An.EO.O. )GU Tfl 90 032 

IMX.G6.0. .AND. AB.6T.10. )GU TO 60 033 

Y A ~ ARS ( Y ) 034 

1 F { X. L6 . 0. .AND. YA.GT. 10. )G0 TO 80 039 

1 F l YA-X.GE. 1 7.9.0R . YA.GE. 6. 9.0R. X+YA.G6 . 9. 9.0R.X.GE. 3. )G0 TO 20 036 

JF(X. 16.-9. ) GO in 40 037 

I F C YA-X.6E.2. 9) GO 1U 90 038 

IF( X4YA.G6. 1.9) GO TO 30 039 

10 N = 6 . *♦ 3 . v AH 040 

615=1. /(N-l. )-Z/N**2 041 

15 N = N - 1 042 

615=1. /<N-1.)-Z*E15/N 043 

1F(N.GF.3) GO Tfl 19 044 

F 15 = 2*6 1 5-CWtXC . 5772 16+ALOG ( AR ) - AT AN? ( Y*X } ) 04 9 

GO TO 90 046 

20 Jl=l 047 

J 2 - 6 048 

GO TO 31 049 

30 J 1 = 7 090 

J2 = 2 1 091 

31 S-I.0,.0) 09? 

YS=Y*Y 093 

DO 3? l=Jl,J? 094 

X 1 = V H J ♦ X 059 

CF =w 1 I )/(X!*X|+YS) 096 

32 S = S*CMPLX(X1 *CF, -ya*CF ) 097 

GO TO 94 098 

40 T 3 = X * X — Y*Y 099 

T4=?.*X«YA 060 

T5=X* T3-YA* ) 4 061 

16=X*14*YA*T3 062 

Fig. 7. Subroutine EXPJ. 
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UC*CMPLX { D( ll)f|.( 12)*X4 0( 13)*T34T6-F( 1?)*YA-E< 1 3 ) * T A , 06 3 

2 EMl)*f : (m*X'fE(i:4)«T3*T64(}n2)*rA + (J(13>*T4l 06 A 

VCsCMPLX(0( iA)4|)( 16 ) *>40(1 to )S=T 34Tt>-H 16)4YA-F< 16)*TA, 066 

2 EU4)+E<16)«X+E( 16)*T34Tf>+l>(15)*YA + l)l 16>*1M 066 

G( l TO 52 067 

50 T3=X*X-Y*Y 068 

T A=?.*X*YA 069 

T6-X< r 7‘ l -YA$TA 070 

T6»X# TA4YA* T 3 07 I 

T 7 = X*T 5-Y A*T6 072 

1 6sX*l 64YA* l 6 073 

T9»X«1 7-YA*TR * OVA 

T ] 0- X* T84 YAt 1 7 076 

UC-CMPLX l 01 1 ) 4 1 ) ( ? ) * X 4 0 ( 3 ) * 7 3 ♦ 0 { 4 ) * T 5 4 1) t 5 ) * T 7 4 T9- ( F ( 2 ) ^ Y A 4 E ( 3 1 * T A 076 

24E( A )<=T64f ( 5 )«» 1 h ) f h( 1 )+E« 2)< K4t ( 3 )< T34h I A JOT646 ( 61 « r /4t 104 07 7 

3 (fH 2 ) * Y A 4|) ( 3 ) “ 1 A4 1) i A ) *> 1 6 4 1) ( 6 ) * I H | ) 07 8 

VC t C M P L X ( f> ( 6 ) 4 1) { / ) *? X-4 0 ( 8 1 * T 3 4 0 ( 9 ) T 6 4 0 ( ] 0 ) T 7 4 T 9 - < E ( 7 ) * : Y A 4 1 ( 8 ) ^ 1 A 079 

2^F ( 9 )*T64E ( 10 )< IF. ) , r ( 6 )4M 7 »^X4E ( b )<‘T34E( 9 )*1 54fc ( 10) ->1 741104 08 0 

3U>( 7)4YA4f)( B)vTA4D<9)*T64|)( 10)*T8) ) 081 

62 FC=UC/VC 082 

S*EC/CMPLX| X, YA ) 083 

6A FXstXP(-X) 08 A 

T-FX*CMPLX(COS( YA),-S1N( YA ) ) 086 

f . 1 6 = S * 1 086 

56 IF(Y.LT.O. )E 15=C0NJG( El 5) 087 

GO TO 90 088 

8 0 F 15 = . A09319/ ( 7 4. 1 930 A A ) 4. A 2 1 83 1 / ( 7.41.02666)4. 1 A7 1 26/ ( 7 4 2. 66788 ) 4 089 

2*2063 35E-1 / ( Z 4A.90C35 )4. 1o 7AC)1E-2/( 7^8. 1821 6 I 4. 1 6865Afc-A/( 24 OVO 

312. 7 3A2 I 4.3 W031 1-7/ ( 74)9. 3967) 091 

t 1S=F16*CEXP( -/ ) 09? 

90 »P(JtN.cQ«l)*)12-Ei5 093 

100 7 =V2 09 a 

7 =V2/V 1 095 

TM = A T AN2 I A I MAG I 7 ) , R E HI 7 ) ) -AT AN2 ( A1 MAO( V2 ) »KEAL ( V2) ) 096 

24 AT AN? ( A I MAC; ( VI ) , RIAL ( VI ) ) 097 

AB=A8SilH) 098 

IFfAB.LT. 1. ) TH ”. 0 099 

1F( TH.GT. 1 . )TH=6. 2831853 100 

1FITH.LT. -1. )TH= — 6. 2R31853 101 

W1?=W12-E154CMPLX( .0, 1H) 102 

RETURN 103 

END 1 OA 


Fig. 7. Subroutine EXPJ - continued 
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SUBROUT INF I DANT ( 1 A, l B* ICC » INS* I NT , I ) « 1 ? • I 3 f J A # JB, JPP, MO, N, NCH, 001 

2ND , Nt 0 , NP , NPGP , NR S . NS • AK , C * CMM , 0 . FMC • CDK , SDK * X , Y , Z • /H , Z l U 1 002 

COMPLEX ZS*ZH»P(2*2)*0(2*2)*CIJ*2SAM*ZUPP 003 

COMPLFX Cl ICC* ICC) *7ll)( 1 ) OOA 

DIMENSION X( 1 ) f YC 1 ), U 1 ), I A( 1 ), 1RU >,NU( 1 ) ,COM 1 ) , SOM lit 0( 1) 005 

DIMENSION I 1 ( 1 ) , I2(!) , 131 1 ) . JA( 1 ) f JB( 1 ) ,MD( INS* A) 006 

DATA TP/6.2B31H/ 007 

2 FORMAT (RX t 'AK= • f FR. 6, 5X, 'DMAX= »,Fd.A f 5X, »DM1N=* ,F8. A) 008 

DO 10 I = 1 » NC M 009 

DO 10 J^l.NCM 010 

10 CI I , J ) = ( . 0* . 0 ) Oil 

0 M A X = , 0 012 

OM! N= 100. 013 

DO 20 J = 1 » NRS 0 1 A 

D J=D( J 1 015 

IF (DJ.GT.DMAX )DMAX=DJ 016 

IF(DJ.LT.DM1N)DMIN=DJ 017 

CDK ( J ) =COS ( f) J ) 01b 

SDK! J ) = SIN( Dj) 019 

K=J*NRS 02 0 

COM K )=CDK( J ) 021 

20 SDK ( K ) =SDK ( J ) 022 . 

I F ( DM 1N.LT.AK ) GO TO 21 023 

IFIDMAX.GT.3. >G0 TO 21 0?A 

IF( AK.GT.O. 1 ) >0 TO 21 025 

GO TO 22 026 

21 WRI TF(6,2)AK t UMAX, DM IN 027 

N~0 028 

RETURN 029 

22 DO 200 K = 1 * NS 030 

NDK—NImK) 031 

K A *= J A I K ) 032 

KH=IB(K) 033 

DK=D(K) 03A 

DO 200 L = 1 * NS 035 

NDL s=ND I L ) 036 

LA& 1 A ( L 1 037 

L f B ( L 1 038 

DL =D( L ) 039 

N 1 L c 0 OAO 

DO 200 I I = 1 1 NDK 0 A 1 

I=MD(K,11) 0A2 

IF( I.GT.NCM) GO TO 200 0A3 

F 1 = 1 . OAA 

I F ( KB. EO. I ? ( I ) ) GO TU 36 0A5 

IFIKB.EO. 1 1(1 ) ) F l = - 1 . 0A6 

1 S= 1 0 A 7 

GO TO AO OAR 

36 1FCKA.E0. 13( I ) )F!=-1. OAV 

I S = 2 050 

AO DO 200 J J = 1 * NDL 051 

J=MD(L,JJ) 052 

1 F ( I . GT. J )G0 TU 200 053 

F J = 1 . 05 A 

!F( LB.EO. I 2 ( J) )G0 TO A6 055 

IFUB.FO. I 1 ( J ) ) F J = - 1 . 056 

JS= 1 057 

GO TO 50 058 

A 6 IHLA.EO. I 3 ( J) )FJ*-l. 059 

J S = 2 060 

50 !F(Nll.NF.O)G() TO 168 06l 

N I L = 1 062 

Fig. 8. Subroutine IDANT. 
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IFfK.EO.l )GI) TO 120 063 

1N0=( LA-KA )* ( LB-KA )< ( LA-K8 J* ( LH-kH) 069 

I F C 1 NO. EG. 0 ) GO TO 80 06b 

C SEGMENTS K AND L. SHARE NO POINTS 066 

CALL 7GSC X( KA ) , V ( KA ) , 7 ( KA ) , X( KH ) , Y( KB ) • Z < KB ) , >. (LA ) , Y ( LA ) , Z ( LA) , 06 7 

2X ( LB ) , Y( LB ) • / ( LH ) , AK, OK, CD*( K ) , SDM K ) « UL t SDK ( L ) . 1 NT , 068 

3P( 1« 1 )«P( 1 ♦ 2 ) * P( 2* 1 ) • P( 2* 2 ) ) 069 

CD TO 168 070 

C SEGMENTS K AND L Shake ONE POINT (THEY INTERSE-Cl) ’ 071 

80 KG=0 072 

JM = KR 073 

JC=KA 07A 

KF = 1 07b 

INO=(KR-LA )*( KR-LH) 076 

IF ( INO.NE.OlGO TO 8? 077 

JC=KR 078 

KF =-l 079 

JM^KA 080 

KG r 3 081 

82 LG=3 ° 8 ? 

JP = L A 08 3 

L F = - 1 08 A 

JFUB.EO. JOGO TO 83 08b 

J P = L B 086 

L F = 1 087 

L G = 0 088 

83 SGN=KF*LF 089 

CPS1 = ( ( X ( JP ) - X ( JC) )MX( JM)-X( JC) )♦( Y( JP)-Y( JC) )*( Y( JM)-Y(JC) ) 090 

?4( z (JP)-Z ( JC ))*(/( JH) -7 ( JC )) )/( DK*DL 1 091 

CALL 7GMM( .o, ok, .o,ol , ak .cdkck ) , SOKI K ) , SOM L ) *CPS 1 09? 

2 , C ( 1, li,w( 1 / , w i ? , 2 > ) 093 

00 98 KK = 1 * ? 099 

KP= I ABS( KK-KG ) 09b 

OU 98 L L = 1 * ? 096 

LP= I ABS( LL-LG ) 097 

P( KP f L P ) =SGN00< KK, LL ) 098 

9 8 CONTINUE 099 

GO TO 16 8 H»0 

C K=L (SELF RFACTION OF SEGMENT K) 101 

1?0 S = . S 10? 

!F(KA. NF. LA»S=-. b 103 

CALL 7GMM( .0,l)K ,0K»( .5-$ ) ,DK»( ,5 + S ) , AK f COK ( K ) , S 0 K ( K ) • SDK (K) , 1 . 1 04 

2, PC 1 , 1 ), *C U ?),p( 2. 1 ), P(2, 2 ) ) 10b 

168 C1.J = F 1*F J*P< I S, JS) 106 

IF ( J.GT.NCM1GI) Tn 190 107 

C( 1 • J ) -C< 1 * J ) 4f. 1 J 1° 8 

!FC l.NE. J ) C ( J« I )=C( J* 1 ) ♦C I J 109 

GO TO 200 1 1° 

1 90 JG = J-JPP 1 1 1 

C( 1* JG)-C( I . JG)-CI J 11? 

200 CONTINUE 1,3 

7M=(.0,.0) 11^ 

)F(CMM. LE.O. )G0 TO 262 H* 

CALL ZSURF(AK f C.MM,FMC,ZS) H 6 

ZH=ZS/( 9. * TP*AK ) 11 ? 

DO 260 K = i ,NS 1 1 H 

NDK=NO( K ) 11* 

ZSAM=2.*/h*(D(K )-SOK(K)*CDK(K) )/SDK(K)**2 1?0 

DO 210 II = 1 , Nl )K 1?1 

I = MD ( K f I I ) 1?? 

210 I F ( I . LE* NCM ICC !• I )*C( I, I ) ♦ZSAM 1?3 

IF (NOK.EO. 1 )G(J 10 260 1?* 


Fig. 8. Subroutine IDANT - continued 
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Zoee-7.., H MS , 1K(K! . 1 , (K) , CUK(K)J/SDMK)oe2 

KA=|A(K) 

1 2b 

KB= 1 B ( K ) 

126 

DO 260 I I a 1 , NDK 

127 

I=MO(K, | | ) 

128 

1H I • GT . NCM )GU 10 260 

129 

FI-1. 

1 30 

JM KB. to. I?( 1 } ,r,o Id p 36 

131 

I F < KB. F 0. 11 ( 1 ) ) E 1 =- J . 

132 

I S- 1 

133 

GO 10 240 

136 

236 IF(KA.60.13( 1 1 )M=-1. 

13$ 

1S = ? 

1 36 

2<*0 DO 260 JJ«l f Nl)K 

137 

J=MO( K, JJ ) 

138 

IF I t . GF. JJGil TO ?60 

1 39 

F J - 1 . 

160 

IF ( KB. F 0. I P( J ) )(;n TO ?46 

161 

IF ( KB. f (J. 11 / J J )f j~-i # 

K? 

JS= 1 

163 

GO 10 2S0 

1U 

2<»6 1 E( K A , E 0 . I 3 ( J ) )FJ»-1. 

14$ 

JS*2 

146 

?!>0 IF( IS.tO. JS)CIJ=FI*Fj*/SAM 

147 

1 Ft !S.NF.JS)CIJ~f- |CF jc?np(> 

148 

IFU.GT.NCM1G0 TO 2b9 

149 

C ( I * J ) -C ( ltJ)*ClJ 

1 $0 

C ( J » I I ( J t I ) ♦c u 

Jbl 

GO TO 260 

1 $2 

2b9 JG= J-JPP 

153 

I • JG 2 -C ( I • J*i> j — C I J 

1 $4 

260 CONTINUE 

lbb 

262 I F ( NL D. 1 E. 0 )G0 TU TOO 

1 $6 

DO 280 I =1 ,NCM 

1 57 

J JA « JA ( 1 ) 

1 $8 

J l ” J J A 

lb9 

I I 2* I 2( I ) 

160 

1 I UI1 ( I ) 

161 

IFf I 12. EO. IHIJ1) J Jl=Jl*rv«S 

162 

IF( 1 . L E . NPQP ) GO 10 270 

163 

JJ8=JR( I ) 

164 

J2*=JJB 

16b 

IF( I 12. FO. Ili( J?) ) J2 = J2*NW<; 

1 66 

C( If 1 )C{ 1 f ] ) ♦ 7 L0( j j ) +1 L(>( J2 ) 

167 

J J J ~ J J A 

168 

DO 268 K c 1 f 2 

169 

NDJ*ND( JJ j ) 

1 70 

DU 266 J J = J , NUJ 

171 

J = Mf)( JJ J, j j f 

1 72 

I F J J. FO. I) GO If) 26 4 

173 

1 F ( I 21 J ), NE. I I 2 ) GO TO 266 

174 

F I r 1 . 

17 b 

I F ( K. EO. 2 )G0 10 2 64 

176 

I F ( I 1 C J ) .Nt . ( I 1 )E I =-J . 

177 

Cl If Jt C< 1 * J 1 + h t * ZH)( Jl ) 

1 78 

GO in 266 

179 

26A 1 F ( I 3 < J ) . NL . I 3 M ) ) E I =~ ) # 

1 80 

C( I* JI«C( 1»J)4F|«ZlD< J2) 

161 

266 CONI r Not 

162 

268 JJJ=JJB 

183 

GO TO 280 

184 

270 IF( 1BCJ1 ) . F .NPGP ) J 1 = J 1 ♦MRS 

18b 

1 «6 


Fig. 8. Subroutine IDANT - continued 
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CM. I >=C( I, I )+2.*2lOi J1 ) 

NOJ=NO( JJA ) 

DO 27 8 JJ = 1,M)J 
J=KO( JJA f JJ ) 

1 F ( J.F.O. 1 1 GO Tf) 2 78 
1 F ( I 2( J ) »NE • 1 I? IC>0 TU 278 
F 1 = 1. 

IF M 1 ( J ) . NC . M 1 IF 1 *-l . 

CC 1. J)=C( 1. J 1+2. *F1*ZLD( J1 ) 

278 CONT 1NUE 
280 CONI INUE 
300 RETURN 
END 

Fig. 8. Subroutine ID ANT - continued 
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192 

193 

.194 

199 

19h 

197 

19 8 

199 
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i , i ? . r 

. SDK , 1 * 
pi,fp2 


MO, N « NO. NS . CDK t C J « 0 , 

X , Y • 7 I 


SllAKflUl INF l PP LIH 1 A, 1 hf 1 NS . i 
26 PP, FT I , Fph, F GPP ,G T i . ph 

COMPLEX I P H » J 1 »«• C.i I » F M » r * 2 « 

i : i » rt r,i i ( , t . ,nK t , , . o, , , . x , , , . v « x 

oi m ns iun m l > . i?< n . i ^ n ) tHOMNS.* > 

f) A I A CJJ/I.Ot-. SiO*8hP-?)/ 

THRs.O) #S*>3 A a 1 H 

C. 1 M*COS ( I I 

srn = s l m c in* ) 

pHRs . 01 7SS *> 4* PH 

cf»H=cns( ph« ) 

SPHsS!N( phk ) 
rx i 1 30 i r i * n 
F 1 T ( 1 ) ■ ( • 0, • O ) 

130 F P P < l ) = ( . 0 ♦ . 0 I 
on iso K= 1 , NS 
KArlMKl 

CftLL^lfMXtKA). V(KA’.?CK A ),*«««». »«*'•« 1 • ? f ^ K ’ 

2 % coki k ) . s(K( k i tC th, s n 


t.ni) 


11, C T2.PP! ,r w ^) 


1 3 1 


200 

260 


001 
00 ? 
003 
oos 
oos 
006 
00 / 
00 « 
00V 
OiO 

01 l 
0)2 
01 3 
01 s 
oi s 
016 
0 1 / 
01 * 
OlV 
020 
021 
02? 


MlKsNfU K 1 


0? 3 

nn iso li-i * nok 


o?s 

1 =HD( K, 1 1 ) 


02S 

M-l. 

1 >6 

026 

1F<KR.F0.1?< 1 1 ™ 

027 

IHKH.frO. UllM* l=-l- 


0?« 

F P P ( 1 1 b F P P C l 1 ♦ I * P p 1 


02 V 

Finu-.mni^^ni 


030 

GO TO ISO 


03) 

1 F 1 K A . t f . • ! : l ! J 1 c 1 ~ - 1 - 


03? 

t pp 1 1 ) =1 PP l 1 > 1 *CP? 


033 

FTT(1)-HW*. 1 ♦ H • € T 2 


O a S 

CONI 1 NUP 


G3S 

FPH=( .0.-01 


036 

F 1 H* (#0**01 


037 

00 ?60 ! - l » N 


03H 

F 1 H»f ? »*♦<. JIIDHIII) 


0 3 9 

6f*M«M'H«r. Jf I 1 fc F p P ( 1 1 


oso 

APPerC Ah > I P PH ) 


OS l 

A 1 1 t-CA' SI F 1H) 


Os 2 

CPP-APP^ A P P t ( 30. • & 1 


0* 3 

GTT*ATT»ATT/< 30* *0) 


OSS 

RfllikN 

6M) 


OSS 

Fiq. 

9. Subroutine iFFLD. 



31 


SUHRHUT I NF 1 SOW I(IA.!H,1CC*ICJ«I NS f l^»<IT c *ll l l? f |3 # ^A, m . 

AX t M | N, v l>, N, Nf M , *:i », NP, \*S • K-* P , 6 S , N W , 6 SOP , DC , > C , VC • / C ) 

(» 1 Mf-NS lHN j SP ( .>()), I>c» l I • / Cl 1 I . Vf. ( 1 | t /C | l ) 

0 I HF MS I l*N I 1 ( 1 ; , I ? ( 1 ),!■}(! ) , j A ' ] | , j H I 1 ) 

DIHFA s ION ! All ) , |H( 1 ) t M)( 1 ) ,FD< INS, A ) 

1 I f )K K A ? ( 9 * , ’ J * . IX, * ! A ( J ) « , 1 > , • M ( J ) * , l } X , 1 * * , 

? 1 / . • ! A ( K J » . 1 X , • 1 * f K ) * , 7 / , * nc I J ) 1 I 

? FlIftKA T c ) K t 319, 1 I 1 9, ?|9, ?M 9. 9) 

3 I (Jkw/.l I lx, J IS, J M 0 . S , l 19, 3 P 1 0 ♦ 9 ) 

A rPRMAHl>,Mb,lllb,6lt») 

S FORMAT ( 1 *«0 ) 

h M’Kma I I | 4X. «>fCl II •, bX, • YC« I I •, NX, «/ C C , 

? 9X. •.!•, AX, *XC| J) • .9X, *YCI J) *,9X, */C< J) • I 
7 FORMA 1 ( 9X , • | • , 3< , * j A • , AX, • JH • , 3X, « | i • t 3X , • | ? • , <U t • I 1 « f 

?1 AX, •* • , 4X, • JA* , iX , • Jji* f 3X, * l 1 1 , .1 X , * | ? 1 , 3 X , • J 3 • i 

1 

NP I = Nwp-».P0.P 

C Ml* XT S* I UP 7 uf l MA (,F SIGMtNTS 
HO 1H Jrl.NKS 
K=j4NWS 

1 A ( K ) S I A ( J ) 

I F ( U( Jl.GT, NPf,P ) | A ( K ) s | A ( J | •* SP | 

IBIKhIHUi 

7 8 1 F ( IHCJI.GT. NMGM ) I H ( K I - I h( J | «NPI 

C NF X 1 SF T l»P Tm£ | MA;,f PUIMS 
no ?0 IsIGPP.NWP 
J * I * NP I 

x c c j ) = x c i n 

YC { J ) * VC ( I ) 

? o zc<j>«-zr.m 

C .\rXT CulCul Air Tmc (. tNf»7«S OCiJl 

IH IWRI TF.lF.O)GU IIJ 
WR IT F (6, 9) 

WR ITH 6, i ) 

?? 00 ?9 J*1.NRS 

K = 1 A ( J | 
t e I h I J ) 

OXrXCC K l-XC ( L ) 

OYbYCC K >-yf ( L ) 
t>Z = /C(K|-/C<L J 

OC I J I -SORT (OX *OX^Oy*()Y4DZ*OZ ) 

Kt J«NRS 

nc I k i -nc ( j i 

?s ! F ( I«^l ! F • OF . 1 )WK | lM6f ?) J, IA(J), |h(JJ,K,|A(KJ,lK(*J,OCIJ) 
I F ( I W« i T f # L 6 • 0 I f*(; T n j? 

WP I TF. U , a j 
WW|TH6,M 

00 30 l=|, NWP 

! F ( I • G T . NRfiP ) (*,() TO pf\ 

WR | n C 3 M , xcc I ) , YC ( I) # ZCI II 
GO T»> *0 

?fl J«hM l 

Wk | T M N, 3 ) 1 ,xcm , YCl I ) , ZCll I • J. XC( J ) , YCI J) , / CC J> 

30 COST l M IF 

WR 1 1 f | 6, S ) 

C CMFCK INPUT OATA FOR CliNSlSTFNff 
3? N«0 

M !N* 100 
MAXs: 100 

IF IMRGP.ie.OIGn TO AO 
no 38 Ici.NPGP 

1 *0 


001 

00 ? 

00 3 
OOA 

oos 

006 

007 

OOA 

009 

010 

01 1 
01 ? 
(•13 
01^ 
«)1S 
016 
01 t 
018 
019 
G?0 
C? I 
0 ?? 
0? * 
0? A 
0?6 
0? 6 
0? 7 
0?R 
0? 9 
030 
U3i 
03? 
03 3 
03a 
039 
036 
03 7 
0 3 m 
(•3v 
OAO 
(A J 
OA? 

DA 3 

OAA 

0A9 

C'Ao 

OA f 
CAW 
GA9 
090 
09 1 
C9? 
09 3 
09 A 
099 
096 
09 7 
09R 
099 
060 
06 1 
06? 


Fig. 10. Subroutine ISORT. 
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00 35 J-ltNSGP 

K=(iAij)-i)*(imj)-n 

35 I T IK. fcO. 0)L=L +1 
/T ( (. . L T . MAX ) MAX = L 
3B N = N-f?«L-l 

A0 I F I NRP. L£ . NPGP )C,fj TO 50 
DO 46 I=1GPP,NKP 
L =0 


IJU 


46 

50 


IT (K.60.0)L=L+1 
IF(L.ll.M|N)M|N=L 
N=N+2* ( l-1) 

IHN. lF.o .OK. W.GT. ICJ)G0 TO 500 


IF <MAX.lt. 0 


.OK . 


I f ( NPOP. i. 6.01 GO TO 58 


M !N.lt.0M>() TO 500 


SET UP THE MODES A I 
DO 56 I=1,NPGP 
J = 0 


1HE GROUM.) PLANE THAI 


S2 J — J 4 1 

Ujs|A( J) 

I BJ - I B( J I 

KK = ( 1 A J - I ) * ( 1 BJ- I ) 

I E( J. EQ.NS'GP IGD TO 54 
IF(KK.NE.O)Gf) TO 5 2 

5* J A ( I ) =J 

JB( 1 ) = J + NR S 

IX 11*1 
1 1 1 I ) - IB J 

IE< I 8 J. EO. I )M(I ) = 1 A J 
^6 I 3 C I ) = I l < I )+NP! 

36 I=NP0P 
N-NPGP 
NCM-NPGP 
JPP = 0 


I E I NRS. EO. NPGP )C;0 TO 75 
C SET UP the re SI ()• THE REAL MODES 
00 f.5 K - l » NR P 
N JK - 0 


DO 60 Jri f NRs 

INI)=( I A < J > -k ) c ( | R ( j )-k ) 

I E ( I NO. Nt . 0 ) GO Td 60 
N JK = N JK 4 1 
JSP( NJK ) - J 
60 CONTINUE 
MOD-NJK- 1 

H IMOD.IF.OIGO TO 65 
DO 62 I M()= 1 f M/JD 

I = I 4- 1 

I P0 = I MD* 1 
JA I - JSP( I MO ) 

J A ( J ) = j A I 
JBI =JSP( IPO) 

JB< 1 ) = JB I 
II ( I ) = I A ( JA 1 ) 

1M I A ( J A I I.EO.K) 11(1 ) = I ft ( J A I ) 

I ? ( I I ~K 

I3( n = I AC JBI ) 

6? I H I A t JBI ) . fc OVK ) I 3 ( I ) = I ft ( JB I 1 
65 CONTINUE 
N’C M = I 


JPP=NCM-NPGP 


Will NOT HAVE 


IMAGES 


063 

064 

065 

066 

067 

068 

069 

070 
07 1 
0 72 

073 

074 

075 

076 

077 

078 

079 

080 
081 
082 

083 

084 

085 

086 

087 

088 
069 

090 

091 

092 

093 

094 

095 

096 

097 

098 

099 
TOO 
101 
102 

103 

104 

105 
2 06 

107 

108 
109 
1 10 
1 1 1 
112 
113 
1 14 
115 
1 16 

117 

118 
1 19 
120 
121 
122 

123 

124 
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c SET UP THE IMAGE MODES 
DO 70 1=1GPP,NCM 
K = 1 ♦ J P P 

J A ( K ) =JA( l ) +Nk S 
JR ; K ) = JM 1 ) +NRS 
1 1 A= I 1 < 1 ) 

110=1211) 

1 I C= 1 3( 1 ) 

11 (K ) = 1I A 

! F C 1 I A . GT.NPGP ) 111 K ) = 1 I A*NP1 
1 2 ( K ) « 1 I H 

1F( 1 IB. GT.NPGP) I 2 C K ) = I 1B + NP1 
1 3 ( K ) = 1 1C 

7 0 1 F ( 1 IC.GT.NPGP ) ! 3(K ) = I 1C + NPI 

N=?*NCM-NPGP 
75 M A X - 0 

M 1 N= 1 00 _ , 

C NO ( J ) = NUMBER OF DIPOLE MOOES SHARING SEGMENT J 
C KD(J,K) = LIST OF DIPOLES SHARING SEGMENT J 
DO 100 J = l » NS 
DO 80 K= 1 * 4 
80 MD ( J » K ) =0 
K = 0 

00 90 I =1 t N 
J A I = J A ( I ) 

JBI=JB< I ) 

L=( JA1-J)*( JBI-J) 

IF(L.NE.O)GO TO 90 
K = K4 1 
M D ( J » K ) = 1 
VO CUNI 1 NUb 
ND( J ) =K 

1F1K.GT.MAX ) M AX =K 
100 IFIK.LT, MINIMI Ns K 

I F 1 1 WR 1TF.LE.0) GO 10 500 
WR I TE ( 6f 7 ) 

DO 1 1 0 I = 1 . NC M 
I F( I . GT.NPGP ) GO TO 1 08 

WR1T616,4)I,JA( 1 ) , J i> (I ) ♦ 1 HI ) * ! 2 ( I ) » 13(1) 

GO TO 110 

1 n8 y - \ + IP P 

WRITE(6*4)1»JA( 1 ) » J B ( 1 )» 11( I )♦ 1 2 C I ) * 1 3 ( I ) • K * J A 1 K ) * J B 1 K ) * 

21 1(K) , I2IK ) , I 3 ( K) 

110 CONTINUE 

WR 1 T E ( 6 * 5 ) 

500 RETURN 
END 


12 5 
126 

127 

128 
129 
1 30 

131 

132 

133 

134 

135 
1 36 

137 

138 

139 

140 

141 

142 

143 

144 

145 
1 46 

147 

148 

149 

150 
161 

152 

153 

154 

155 

156 

157 
168 
159 
] 60 
161 
162 

163 

164 

165 

166 

167 

168 

169 

170 

171 


Fi.g. 10. Subroutine ISORT - continued 


34 


/ 


SUBROUT INF P 0 ! S S ( 1 A , IB, INS, I i « 12. I 3, MO, NO, K'LD, NS , C J , CMM, o, COK , 
2SDK,OISS,ZH,ZL')l 

COMPLEX ;■ J( 1 ) t ZH t C.IA, CJB, 2LA, 2LB. 2CO( 1 ) 

DIMENSION CUKIl ),SUK( 1 ) , 0 ID , 111 I ) , 1 2 ( 1 ) , I 3 ( 1 ) , I A ( 1 ) , 1 B ( 1) , NO ( ] ) 
DIMENSION M 0 ( INS, 4 ) 

RN=RF AL ( 2H ) 

DlSS=.0 
00 100 K= 1 ,NS 
IFICMM.LE.O. )GU TO 10 

EA = 2.*KH*(()(K ) -SDK ( K )*COK( K ) ) /SDK! K 1**2 
F B = 4 . *RN* I SDK ( K ) -01 K ) OK I K ) ) / SDK I K ) * * 2 

1 0 K A = 1 A ( K ) 

K B = I B ( K ) 

C J A = ( .0, .0) 

C JB=( . 0, . 0 ) 

NOK =N0 ( K ) 

DO AO 1 1=1, NOK 
I = MO ( K , i I ) 

F I = 1 . 

IFfKB.EO. 1 ?{ I ) ) GO TO 36 
1FIKB.EQ. I] I I ) )E I = -l. 

CJA=CJA+FI*CJ( I ) 

GO TO 40 

36 IFIKA.FO. I 3 ( I ) )F!=-1. 

CJB=CJB«F I*CJ< I ) 

40 CONTINUE 

1 1 ( NLD.LF.OIGO TO 60 
A JA=C AhS ( C JA ) *<■ 2 
BJB=CABS( CJB )**2 
K K = K ♦ N S 

RIA=REAL I 2LDIK ) I 
RLB = RFAl (Zll)( KK ) ) 

0 I SS = D I SS+A JA*KL A + HJB*RLB 
60 1 F ( C MM . L t . 0 * ) GO 10 100 

DISS=DISS+Fa*<CABS(CJA)**?*CABS(CJB)**2) 

2*FB*I RE AL ( C JA )*REAL ( C JB 1 +A 1MAGI C J A )* A I MAGI C JB ) ) 

1 00 CONT 1 NUF 
RETURN 
END 


001 

00? 

003 

004 
006 
006 
00 7 
008 

009 

010 

011 

012 

013 

014 
016 
016 
017 
01P 

019 

020 
021 
0°2 

023 

024 
026 
026 

027 

028 

029 

030 

031 

032 

033 

034 
036 

036 

037 
03B 
039 


Fig. 11. Subroutine PDISS. 
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SUBROUTINE ZEE ( XA, YA, ZA,XR*YB,Zft*U 001 

2tCK0, SKOt CTH f S rn, CRH, SPH, f T 1 , F T 2 , EP 1 • EP2 ) 00? 

COMPLEX EJA,EJB,EPUEP2tESl ,PS2tETl f ET2 003 

C A = ( XB-XA )/l> * 00^ 

CB-(YB-YA)/D 00b 

CG= ( Z B-2 A ) /D 006 

G=(CA*CPH+CB* SPH ) STH+CG^C TH 00 7 

GK=1.-G*G 008 

ETl=(.Ot.O) 00V 

E T ?= ( • 0» • 0 ) 010 

EP1=(.0,.0) 011 

f:P2-(.O,.01 012 

I P t GK ♦ L T • • 00 1 ) GO TO 200 0 * 3 

A = XA* si H*f.PH + YA* STM® SPH4ZA*CTH 01 A 

B=XB*S1 H*CPH4>YH*S * M *SPH*Z H*CT.H 01b 

E JA=CMPLX(CnS( A ) , S 1 N ( A ) ) 016 

E J B = C M P L X ( C 0 S ( B ) , S 1 N ( B ) ) 017 

SGD=SIN(G*D) 018 

CGD=COS(G*D) 019 

ESI ^30.»EJA^CMPl X{ SGO-G»SKO,CKD'CGO)/GK/SKO 020 

ES2=30.*E JB^CMPLX ( G^SKD-SGO? CKO-CGU ) /GK/SKD 0? 1 

T= ( CA*CPH+CH*SPH ) *CTH-CG*STH 022 

P=-CA*SPH+CB*CPH 023 

ET1=7*ES1 02A 

ET2=T*ES2 025 

EP1=P*FS1 026 

FP?^P*ES 2 027 

200 RETURN 02ft 

END 029 


Fig. 12. Subroutine ZFF. 
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a 


3 




Til Ti 2, T01 1 102* SO! tO^SItODt ZO 

COMPLEX FGZ»2,2J.r.M<2».GIM2> 

DA I A F T A, GAM, P 1 / 376. 727, t . 0. 1 • 1 • 3* 1 1 - ,q - 

DSO=D*D 

SGDS=SGOl 

1 F ( S2 % l T .SI ) SG0S-~SGD1 


SCO T = SG02 

|F(12.Ll.Tl) SGO 1 =- SGD? 
IF( AhS(CPS!).GT.O. 997)00 
ESl=CFXP(GAM*Sl ) 

E S2 = CF X P ( GAM* S2 ) 

El 1=CFXP(GAM*1 1 ) 

ET2=CEXP(GAMoT?) 

00 = 0 

0PS1=CPSI 
TO 1 =T I 
TD2 = T 2 

CPSS=DPS1*UPS1 

CD=00/0S0RT ( 1.00-CPSS) 

C = CO 

B0=C0*0PS 1 


T(1 


110 


E B = CE XP ( G AM^CMPL X ( . 0 , B ) ) 

EC = CE XP l GAM*CHPLXC • 0, C ) ) 
00 10 K= l » 2 
00 10 L = 1 * 2 

10 = 

tsi=toi*toi 
TS 2=T02*TD2 
OPO = DO s> 00 
S I =S 1 

00 100 1 = 1 f 2 

SOI=SI 
S 1 S = SD 1 * SD I 
ST1«2.*SD1*TD1*0PSI 
ST2 = 2. <r SOI*TD2 <s OPSI 

R 1 =DSQK 1 ( DPO**- S l S+T S 1-ST 1 ) 
R2=0S0RT ( DPO + S IS+1S2-ST2) 


40 

SO 


= EH 

SO K = l » 2 
= ( -1 )**K 
= F K* SO 1 
= EC 

4,0 1 = 1.2 
= 1-1 ) * *L 
l=FK*Fl 
sFK^BO+FL^CD 

l=R_*TDl 

2=Fl*T02 

tl =R USK + Tl 1 

[K.D-tlK.Ll + H* ( EXPA^EKL + FXPB/EKL ) 

_= 1 . /FC 
K -\. /E B 
n=SPl*DPSl 
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Fig. 13. Subroutine ZGMM. 
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2 C = Z D 063 

FGZ I =CEXP{ GAM*ZC ) 064 

RR1=R1 + ZI)-TU1 066 

RR?=R2-fZ()-TD2 066 

CALL EXPJ ( G4M*RR l,GAM*KR2,EXPB) 067 

RR l =R 1-ZD+T 01 068 

RR2 = R2-ZI) + T02 069 

CALL FXPJ( GAM*Rkl,GAN.*RR2. tXPA ) 070 

F< 1 , 1 ) = 2. *SGDS*( . 0. 1 . )* EXPA/EGZ 1 07 ) 

F ( I, 2)=?.* SODS* (.0.1. )*EXPI1#EGZJ 07? 

100 Sl=S2 073 

C S T =-E T A / ( 16. #PI*SGDS*SGDT ) 074 

P 1 1 =CST # ( ( F( 1, 1 )+E(2t2J*8S2-E( 1,2)/FS2)*FT2 07 5 

A +(-F( 1.2I-H2, 1 )-ES? + E( 1, 1 )/ES2)/M2) 076 

P12=CST#( (-F( 1, l )-E( 2, 2) *FS2+E< 1, 2) /ES2)*E 71 077 

B ♦( F( 1 , ?)+F(2. 1 >*ES2-E( 1* 11/FS2 1/ET1 ) 078 

P?1=CS7«( f~F{2, 1 )-H2,?)*E51+Ff 1,7) /ESI )*ET2 079 

C ♦( F( 2» 2 ) + 6 ( 2 » 1 ) # ES 1-E ( It 1 )/ESl 1/ET2) 080 

P22=CS1*(( F(2, 1 H-E(2,2)*ES1-E( 1.21/ES1 )*ET1 081 

D 4(-F(? t 2)-h(2,l)*ESl+E(l,l)/F$l)/ETl) 082 

RETURN 063 

110 1FCCPS1.LT. 0. )G0 TO 120 084 

T A = T 1 085 

T 6 = T 2 086 

GO TO 130 067 

120 1 A = -T 1 080 

TB = -T 2 089 

SGOT =-SGOT 090 

130 SI=S1 091 

DO 150 1=1,2 09? 

T J = T A 093 

DO 140 J = 1 .2 094 

Z I J = T J-S l 096 

RESORT (DSQ+ZIJ*ZIJ ) 09 6 

W=R+ZIJ 097 

1F(Z1J.LT.0. ) W=OSO / ( R-Z I J ) 098 

V=R-7 I J 099 

1F( Z1 J.GT.O. ) V = OSO / ( R+Z I J ) 100 

i f ( j.ro. i ivi=v ioi 

1 F C J.EO. 1 ) W 1 =W 102 

EGZ ( f t J )-CEXPCGAM*Z I J ) 103 

140 T J = T B 1° A 

CALL EXPJ ( GAM* VI • GAM#V, GP( I 1 ) 105 

CALL EXPJ(GAM*W1,GAM*W,GM< I )) 106 

150 S 1 =S2 107 

CST=ETA/(R.*Pi*SGDS*SGUT ) 108 

PinCST#(GM( 2)-hGZ(?,2)4GP(2)/EGZ(2,2) 109 

2-CGDS v ( GM( l )*fcGZ( 1,21+GPC 1 >/FGZ( 1,2) ) ) 110 

P 1 2 = C S 1 * ( -GM(2)*CGZ ( ?» 1 )-GP< ? l/EGZ ( 2, 1 ) HI 

24CGDS* ( 6M( 1 )*EGZ( 1 , 1 )+GP( 1 l/EGZ (1,1)1) 112 

P2 1 =CST* C GM Cl) *FGZ< 1 , ?) 4GP( 1)/FGZ( 1 , 2) 113 

2-CGDS* (GM( 2 )*EGZ( 2,2>+GP( ? l/EGZ C 2, ? 1) ) 1 1 

P?2 ~C ST* ( -GM ( 1 )*EGZ( 1, 1 ) -GP ( 1)/FGZ(1,1) 115 

2 +CGDS* ( 2 )*EGZ( 2, 1 )4GP( 2 ) /EGZ (2,1))) 116 

RETURN H ? 

END 118 

Fig. 13. Subroutine ZGMM - continued 
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SUBROUTINE ZGSCXA,YA,ZA,X8.YB,ZR.Xl,Yl,Zl.X2,Y2,Z2,AK, 

OATfl L ElA^GAMi^’/376.727,(.0,l. >,3.1M59/ 

C A = ( X2-X i ) /OT 
CK=( Y2-YI ) /OT 
CG= 17?“/ 1 >/OT 
C AS= ( X R-X A ) / DS 
CBS=( YB-YA ) /OS 
CGS=( ZB-ZAI/OS 
CC=CA*CAS+CB*CBS*CG*CGS 
IF tABSICCI.GT. 0.997)00 TO 200 

?0 SZ=1X1-XA)*CAS*(Y1-VAI*CBS+IZ1-ZA)*CGS 

1 F ( INT.tO.OlGU TO 300 
CGOS=COS 

SGOS = Ct*BLX( . 0» SOS ) 

SGOT sCHPl.X I . 0 f SOT ) 

INS=?* ( INI/?) 

1 F ( 1NS.L T . 2 ) 1 NS =2 
1 P = INS* 1 
delt=ot /INS 
T = . 0 

DSZ=CC*DELT 

P 1 1 = ( . 0. . 0 ) 

P12=( .0, .0) 

P2l=( . Of . 0) 

P2? = ( .0, .0) 
aks=ak* ak 

SGN = - 1 . 

DO 100 I N«lt IP 
Z?1=SZ 
Z Z2 = SZ-0S 

XXZ=X1+T«CA-XA-SZ-CAS 

YYZsY1+T*CB-YA-SZ*CBS 

Z2Z=Z1*T*CG-ZA-SZ*CGS 

RS=XXZ**2*YYZ**2+ZZZ**2 

R 1 = SORT ( RS + Z Z 1**2 ) 

EJA=CMPLX(C0S(K1 1.-SINCRIH 

EJl=6JA/P.l 

R2 = S0RT (RS+ZZ2**? ) 

E JB=CMPLX » COS ( R2) * -SI NIR2 ) ) 
Ii(?:wSisOOS*ZZ|.EJl*CGDS-ZZ2.6J2 

eR2=-EJB*SGOS+7Z2*tJ2*CGOS-721-EJl 

lFIRS?GUAKSIFAC-CCA*XXZ4CB*VYZfCG*ZZZI/RS 

ETl-CC*(fcJ2-tJl*CGDS)*Ft *ER1 

6T2=CC* ( E Jl-6J2*CGt»S )+FAC»ER2 
f=3. +SGN 

lftIN.tO.I .UR. IN.EO. IP)C=1. 

C1=C*S1N10T-T ) 

C?=C°S I N( T I 
PI 1- PI I ♦El l* c 1 
PI 2=P1 2+ET 1*C2 
P? 1 =P? i*FT ?*c 1 
P?2=P22*fcT 2*C2 
Tr T ♦OELT 
sz=sz+osz 

1°0 SGN=-Sf^ ^ ,* e iAaObLT/( 12.*P I •SGOS*SGOT ) 

PI 1=CST*P1 1 
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Fig. 14. Subroutine ZGS. 
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P12=CST«P12 
P21=CST*P?1 
P22=CST*P 22 
RETURN 

200 SI 1 = ( X 1-XA )*CAS* ( Y1 -YA ) *CBS«-( Z 1 -Z A }*CGS 

A -S ^ l*CASI*»2-»-fYl-YA-silifCBS)**2 + (Zl'-ZA-SZi *CG$ ) »*? J 
0»K.= ^RHll < wT^?‘ SZ^ * CAS, ** 2 + <Y2 ' YA ” S2^ * CBS, ** 2+<2^ ~ 2A " SZ2 * CGS, '* ^, 

IF(D(lK.Gr.?0.*4K . AND. INT.GT.OIGO fO 20 

I F ( OOK. L T. AK ) ()0K = AK 

RF TukN^ MM * • °* » SZ 2, DDK , CDS » SUS » SDT , 1 . , P ) l*P12*P21 , P?2 ) 

300 SS=SORT< ] . -CC*(.C ) 

CA0 = ( Cr,S=>CB-CrtS*CG)/SS 
chu= ( c As*cn-cc»s*c a ) / ss 
CGO= ( CBS^CA-CAS^CB ) / ss 

DK=( Xf-XA)*CAD+( Yl-YA )*CBl)4( Z 1-ZA)*CGD 
0K=- AHS ( OK ) 

IF (OK. L T . AK )OK = AK 
XZ-XA+SZ*CAS 
YZ*YA+SZ*CBS 
ZZ=2A*$Z*CGS 

XP1 =X 1 -DK*CAD 
YP 1 = Y 1 -DK**C 80 
Z P 1 = Z 1 -DK$CGi) 

CAP=CBS*CGD-CGS*CRD 

CBP=CGS*CAD-CAS*CGU 

CGP=CAS*CQl>-CBS*CAO 

P1=CAPc(xPI-XZ) 4C8P*(YP1-YZ ) +CGP* ( Z P 1-Z Z I 
Tl-rj/SS 
S 1 =T 1 *CC -SZ 

CALI ZC,MM(S1,S1 + 0S, T 1, H+OT,UK,COS, SOS, SOT , CC , R 11 , R 1 2 , P2 1 , P22) 

K t I I J K N 

END 
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Fig. 14. Subroutine ZGS - continued 
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.zs> 


50 


100 


MJPRnUT I \*t- JSUKf ( AK # CMM f FMC 
CHMPIFX BFS.HF51./S 
OAT A FT A * SOT . TP/ ' 7/, v •> >•>-/ , 

JF ( X. f,T. H # )(;n Til 5c> 

T ~ X / P , 

T ? - T ~ 7 
TA=T?*T? 

- 

•■ ;. ' ; : ; ?! ;• ?» > .• ••• •• ; ' 1 ’ ‘ I4_ ws - j "-. wowu, , 

*'j: »» • V " ' 7 * 5 : » ’* W . 31 ) 6751 ,* 

.«► L X f fifcp, fit II • ' 

or?, 1 ,? ,'oo 1 07 * c "“ t> ' Hf Rp - 8t 1 Pf^+fU l P , 

X R - . 70710 fiHlfcx 

XU)./> 

T = l I -.0660355 ■- <l-.0675?Ilf- nil J! * 8 ' iS,CX1+1 - 
BESrf-eCMPt x r Ci IS ( 1 , , s,,.;, , , * yHiuft 35 J*X1-. 3926 v907*xP 

Bt SI i- OCM1-L X , C0S m , si\,n „ ' >0i> 7809 72, X P 

RF1 URN^ ***' ,- ^'* t *" 7 A*BE S/BES 1 /SOT/SOSWE 


ri S- 15. Subroutine ZSURF. 
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